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ABSTRACT 

We have begun a survey of the chemical and dynamical properties of the Milky Way disk as traced by open 
star clusters. In this first contribution, the general goals of our survey are outlined and the strengths and limita- 
tions of using star clusters as a Galactic disk tracer sample are discussed. We also present medium resolution 
{R ~ 15,0000) spectroscopy of open cluster stars obtained with the Hydra multi-object spectrographs on the 
Cerro Tololo Inter- American Observatory 4-m and WIYN 3.5-m telescopes. Here we use these data to deter- 
mine the radial velocities of 3436 stars in the fields of open clusters within about 3 kpc, with specific attention 
to stars having proper motions in the Tycho-2 catalog. Additional radial velocity members (without Tycho-2 
proper motions) that can be used for future studies of these clusters were also identified. The radial velocities, 
proper motions, and the angular distance of the stars from cluster center are used to derive cluster member- 
ship probabilities for stars in each cluster field using a non-parametric approach, and the cluster members 
so-identified are used, in turn, to derive the reliable bulk three-dimensional motion for 66 of 7 1 targeted open 
clusters. The high probability cluster members that we identify help to clarify the color-magnitude sequences 
for many of the clusters, and are prime targets for future echelle resolution spectroscopy as well as astrometric 
study with the Space Interferometry Mission (SIM Planetquest). 

Subject headings: Galaxy: open clusters and associations - Galaxy: fundamental parameters - Galaxy: Struc- 
ture - Galaxy: Dynamics 



1. INTRODUCTION 

1.1. Galactic Kinematics Using Open Clusters 

Open star clusters have long been exploited as tools for un- 
derstanding Galactic interstellar dust (e.g., Trumpler 1930a,b; 
Clayton & Fitzpatrick 1987; Duti-a & Bica 2000), the age 
of the Galactic disk (e.g., Janes & Adler 1982; Twarog & 
Anthony-Twarog 1989; Phelps, Janes, & Montgomery 1994; 
Phelps 1997; Chaboyer, Green, & Liebert 1999; CaiTaro 
1999), the Galactic disk metallicity distribution and age- 
metallicity relation (e.g., Twarog 1980; Friel & Janes 1993; 
Friel 1995; Twarog, Ashman, & Anthony-Twarog 1997), and 
of course stellar evolution (e.g., Sandage 1957; Cannon 1970; 
Maeder & Mermilliod 1981; Meynet, Mermilliod, & Maeder 
1993; Koester & Reimers 1996; Prada Moroni & Straniero 
2002). The value of open clusters as tracers of the local Galac- 
tic rotation curve has also long been recognized (e.g., Hron 
1987; Scott, Friel, & Janes 1995; Glushkova et al. 1998; Lok- 
tin & Beshenov 2003; Frinchaboy 2006a). It is in this role as 
a dynamical tracer of the Galactic disk that the present study 
of open clusters is especially focused. 

Star clusters can be effective tracers of the Galactic disk 
because they offer many advantages over other tracer candi- 
dates. First, relative to other tracers, star clusters lend them- 
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selves more amenably to age, metallicity, distance, and ve- 
locity evaluation. Compared to an isolated field star at the 
same location in the Galaxy, these quantities are much eas- 
ier to establish in a star cluster, needing only a properly in- 
terpreted color-magnitude diagram to establish the first three, 
while velocities can also be better determined for a star clus- 
ter because: (1) averaging radial velocity and proper mo- 
tion data over an ensemble of co-moving stars confers po- 
tentially as much as a ^/N increase in precision for the bulk 
motion of the ensemble, and (2) better distances allow one to 
translate proper motions into transverse velocities more ac- 
curately. The supplemental knowledge of age and metallic- 
ity of a source confers additional beneficial insights into its 
proper use as a dynamical tracer with respect to, for example, 
assumptions about orbit shape and asymmetric drift. Alter- 
natively, one can explore Galactic dynamics as a function of 
population age and metallicity if all relevant data are avail- 
able. 

On the other hand, there are some complications in the 
use of open clusters as dynamical tracers. The challenge of 
proper identification of cluster members can present particu- 
lar hazards. For example, Frinchaboy (2006b) showed that 
the UCAC stars used by Dias et al. (2006) to estabUsh the 
proper motions of at least two particular clusters — Be29 and 
BH176 — in their exhaustive survey of over 400 systems are 
too bright and cannot be part of these distant systems. Dif- 
ficulties caused by inaccurate membership censuses are why 
continued large-scale observational efforts are needed before 
we can be confident in the use of open clusters as Galactic disk 
tracers, particularly for more sparse and more distant systems. 

To overcome these types of problems, which are typically 
associated with small number statistics, it is desirable to sur- 
vey large numbers of potential open cluster members. How- 
ever, this desire to achieve the largest possible statistical sam- 
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pies often encourages a dangerous reliance on compilations 
of disparate data. For example the dynamical study of the 
Galactic disk using open clusters by (Hron 1987) found that 
in their compilation of data from the literature ^ 50% of the 
clusters with multiple distance measurements had differences 
greater than one magnitude in determined distance modulus 
and ^ 50% of the clusters also had poor RV qualities.^ 

To help overcome both the membership and homogeneity 
problems that are often a hindrance to the use of open clus- 
ters as tracers of the Galactic disk, we present a new survey of 
open clusters that will not only take advantage of quahty ra- 
dial velocities to help discriminate cluster members, but also 
rely on one source of data as much as possible for each in- 
dependent cluster parameter (e.g., all photometry from one 
source, aU RVs and derived in a uniform manner, all proper 
motions coming from one catalog, etc.). A central objective of 
this study is the use of these clusters to derive global dynam- 
ical properties of the Milky Way disk, with a particular em- 
phasis on the derivation of the full space velocities for the tar- 
get clusters. In keeping with our philosophy of uniformity of 
data, and because large numbers of proper motions presently 
must come from all-sky astrometric surveys, we first inves- 
tigate the main proper motion catalogs available for investi- 
gations of open cluster kinematics and from which we might 
draw an initial target sample. 

1.2. AU Sky Proper Motion Surveys: Hipparcos, Tycho-2, 
UCAC-2 and the 4M 

A key advance that has propelled a resurgence in the use 
of open clusters as disk dynamical tracers is the compila- 
tion of all-sky proper motion surveys. There have been four 
surveys to determine bulk Galactic cluster kinematics by av- 
eraging the proper motions of presumed cluster stars based 
on data from the Hipparcos (Baumgardt, Dettbarn, & Wielen 
2000), Tycho-2 (Bias, Lepine, & Alessi 2001, 2002a), Four 
Million Star (4M) (Glushkova et al. 1996), and recently the 
UCAC-2 (Dias et al. 2006) catalogs. The ESA Hipparcos mis- 
sion has provided critical astrometry for two of these proper 
motion databases via each of Hipparcos' primary data prod- 
ucts: (1) the Hipparcos catalog of ~ 118,000 stars (V < 11) 
with proper motion uncertainties of 1-2 mas yr"^ and (2) the 
Tycho-2 catalog of 2.5 million stars (V < 13.5) with proper 
motion uncertainties of 1-3 mas yr"' . The UCAC-2 catalog of 
48 million stars is based on Tycho-2 and fainter, ground based 
observations (to R=\6) and has proper motion uncertainties 
of 1-7 mas yr"' . The 4M catalog (Volchkov, Kuzmin, & Nes- 
terov 1992) was compiled from the Astrographic Catalog and 
the Hubble Space Telescope guide star catalog (GSC) reduced 
to the older system of the PPM (Roser & Bastian 1991) sur- 
vey, with proper motion uncertainties of ^ 10 mas yr"'. 

Proper motions for hundreds of open clusters have been 
published using these catalogs. However, a comparison of 
the derived motions for the clusters in common between sur- 
veys reveals substantial discrepancies, as shown in Figure 1 , 
where a cluster by cluster comparison of the proper motion 
differences is illustrated, namely for (a) Tycho-2 versus Hip- 
parcos, (b) Tycho-2 versus 4M, (c) 4M versus Hipparcos, and 
(d) Tycho-2 versus UCAC-2. As may be seen, differences in 
derived proper motions typically exceed the quoted uncertain- 
ties claimed by each survey. The best correlation of derived 

** It is worth pointing out that this is a common problem with other types 
of Galactic rotation tracers that have been adopted in the past, not just open 
clusters. 



proper motions is between surveys using the Tycho-2 and 
Hipparcos astrometry; remaining differences in derived mean 
cluster proper motions between these surveys must therefore 
be due to differences in the adopted samples of presumed 
cluster members because the actual proper motions, at least 
for y < 11, are the same (i.e. HIPPARCOS-based), while 
the Tycho-2 astrometry used at fainter magnitudes is on the 
Hipparcos reference system (which means that the system is 
referenced to background, extra-galactic sources of the In- 
ternational Celestial Reference System). Of course, with its 
bright magnitude limit, Hipparcos can usually provide useful 
astrometry for only a small number of stars per cluster (typ- 
ically less than four). With only a few stars per cluster, a 
Hipparcos-based survey is far more susceptible to small num- 
ber statistics as well as the misidentification of true cluster 
members against the large number of fore/background stars 
of the Galactic disk. 

1.3. A Closer Look at the UCAC-2 and 4M Catalogs 

Clearly Tycho-2 and Hipparcos, which are currently the 
most accurate all-sky astrometric surveys, must be considered 
primary and important sources of proper motion data for our 
survey. On the other hand, deeper catalogs can provide more 
cluster members, but typically with worse precision. Thus, it 
is not immediately obvious that adding additional data from 
the deeper proper motion catalogs improves or degrades those 
from Tycho-2 and Hipparcos alone. A reasonable correlation 
of derived cluster motions is found when either the Tycho-2 
and UCAC-2 catalog data are used with a V < 13 limit, but 
this is because UCAC-2 adopts Tycho-2 proper motions for 
stars brighter than about V = 13 (Zacharias et al. 2004). On 
the other hand, it is clear that there are greater deviations in 
derived proper motions when we incorporate the fainter stars 
from UCAC-2. Apart from not knowing whether these differ- 
ences reflect systematic problems in the fainter UCAC catalog 
(which allows probing of proper motions with stars to V = 16) 
or small number statistics in the brighter surveys, some addi- 
tional concerns about UCAC beyond those suggested by Fig- 
ure 1 have led us to not adopt this dataset for our own work. 

For example, as recently pointed out in Balaguer-Ntinez, 
Galadi-Enriquez & Jordi (2007), the UCAC-2 proper motions 
may have systematic trends with magnitude due to the com- 
piled nature of the UCAC-2 survey (i.e., ground based proper 
motions are added to Tycho-2 data). This concern is use- 
fully illustrated by looking at the cluster M67 (NGC 2682). 
In Figure 2a, we show the 2MASS color-magnitude diagram 
(CMD) for the M67 field plotting only the most probable 
members based on CMD location. We split the CMD (red: 
UCAC mag < 13.0, blue: UCAC mag > 13.0) at a magnitude 
that represents approximately the transition within the cata- 
log from Tycho-2 to ground-based observations. Substantial 
proper motion shifts are apparent between the bright and faint 
samples (Figure 2b), and this suggests significant systematic 
zero-point offsets within the UCAC-2 database. 

The 4M catalog, as well, appears to have systematic proper 
motion errors. The 4M is not tied to the Hipparcos system, but 
rather to the PPM (Gulyaev & Nesterov 1992). Glushkova 
et al. (1996, 1998) have used the 4M catalog to determine 
the proper motions of about 200 open clusters. Dias et al. 
(2001) compared their Tycho-2 open cluster proper motions 
to both those based on Hipparcos (Baumgardt et al. 2000) and 
the Glushkova et al. 4M work and found that the 4M mo- 
tions were systematically offset from those in the Hipparcos 
system by 5 mas yr"' in both p,aCos5 and p,s, an amount 
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that was larger than expected given the quoted errors of both 
the surveys. While these differences likely reflect both dif- 
ferences in membership as well as astrometric accuracy, this 
comparison suggests that the deeper proper motions are not 
necessarily providing better overall accuracy in the open clus- 
ter bulk motions, and recommends a strategy based on quality 
over quantity of cluster star motions. 

Therefore, because of uncertainty over the reliability of the 
UCAC-2 and 4M surveys and our desire to adhere to a "qual- 
ity over quantity" poUcy, we have elected to focus on deriv- 
ing bulk motions using astrometry from the Tycho-2 catalog, 
but with dedication to ensuring that we derive a Irustworthy 
membership of the smaller number of available cluster stars 
available in this shallower database. 

1.4. A New Galactic Tracer Survey 

The mass and mass distribution of the Galactic disk has 
been a matter of debate for over a century, and will hkely 
remain so until extremely precise proper motions and trigono- 
metric parallaxes can be obtained for numerous disk trac- 
ers, most likely through future space-based studies like the 
National Aeronautics and Space Administration's (NASA) 
Space Interferometry Mission (SIM PlanetQuest) and the Eu- 
ropean Space Agency's (ESA) Gaia satellites. The new 
project presented here represents both a preparatory effort in 
this space-based direction as well as a standalone dynami- 
cal study in its own rite. Our goal is to establish a well- 
constructed, well-studied, baseline tracer population — open 
clusters — that can not only (1) serve as input targets for 
Galactic dynamics studies with SIM PlanetQuest (specifi- 
cally, for the SIM Key Project Taking Measure of the Milky 
Way, for which SRM is the Principal Investigator and which 
has provided support for this project), but which (2) can also 
be immediately exploited for understanding Galactic dynam- 
ics with existing astrometric data. 

As mentioned above, the inability to establish a uniform, 
unbiased tracer sample has been one of the key weaknesses 
of previous Galactic dynamical surveys (e.g., Fich, Blitz, & 
Stark 1989). To provide a homogeneous set of tracers, we 
have undertaken a spectroscopic survey to obtain precision 
RVs of open cluster fields. These RVs will establish cluster 
membership for individual stars that not only provides a very 
precise mean RV of each cluster, but, in identifying cluster 
members having accurate astrometry, can be used to define 
the bulk cluster proper motion. The combination of the newly 
found, very precise mean RV of each cluster with its derived 
bulk proper motion and distance will allow us to determine 
the space velocities of these clusters. With a large number of 
cluster space velocities, the rotation curve of the Galactic disk 
can be constrained over the Rgc range of the sample. Alterna- 
tively, through the adoption of an assumed rotation curve (i.e., 
Galactic potential), the orbital properties of individual clusters 
can be determined. 

Because we are interested in obtaining results before SIM 
PlanetQuest and Gaia are in service, our RV study will fo- 
cus on clusters already having available, uniform and reliable 
proper motions. As described in §1.4, we have elected to fo- 
cus on the all-sky Tycho-2 proper motion catalog, which pro- 
vides useful astrometry for typically 50-200 stars per cluster 
field (< 0.75 deg^). While selection of Tycho-2 as our source 
of proper motions limits the depth and thereby the cluster dis- 
tance that can be explored, it is in keeping with our philos- 
ophy of quality over quantity for the astrometric data. The 
selected proper motion stars for a given cluster field can usu- 



ally be investigated spectroscopically with a single pointing 
of the NOAO Hydra multi-fiber spectrographs on the CTIO 
4-meter and WIYN^ 3.5-meter telescopes, and the radial ve- 
locities derived from these spectroscopic data are the primary 
results presented here. Our campaign of multi-fiber spec- 
troscopy allows us to check virtually every star in a cluster 
field having a Tycho-2 proper motion, and leaves additional 
fibers to (1) expand the RV membership census to fainter stars 
in anticipation of the future astrometric surveys (e.g., SIM and 
GAIA), and (2) improve age-dating of the clusters through 
CMD-isochrone fitting to established member stars. The cur- 
rent study of clusters provides a large uniform database for 
further open cluster research, as a supplement to the Dias et 
al. (2002b) and WEBDA (MermiUiod 1995) databases. 

The new RVs immediately improve all previous proper mo- 
tion work on our targeted clusters because of the clarity they 
bring regarding cluster membership. The improved RVs and 
proper motions, when combined with new distances we shall 
derive elsewhere (Paper II), will provide much more reliable 
space motions of numerous open clusters over a large Rgc 
range; these space velocities will be at a precision sufficient 
to make tangible improvements in the determination of the 
nearby Galactic rotation curve and, in turn, the mass distribu- 
tion of the Galactic disk. With uncertainties of order ~ 1 .2 
km s~\ the data here yield the best derived bulk RVs thus far 
for most of the chosen clusters. This precision is comparable 
to the uncertainties in transverse velocity that SIM and Gaia 
will measure for these clusters, and represent a significant im- 
provement over many previous RV surveys of open clusters, 
which have typical uncertainties of order ~ 15 km s~' (Scott 
et al. 1995). Our results are more comparable to the RV pre- 
cisions being obtained for open cluster stars in studies using 
CORAVEL (e.g., Mermilliod & Mayor 1989, 1990), for ex- 
ample. 

Following the work in this contribution (Paper I), we will 
provide uniformly-determined distances and ages derived 
from isochrone-fitting to 2MASS photometry of these clus- 
ters, aided by the cluster membership data derived here (Pa- 
per II). With newly-derived kinematics and distances in hand 
from Papers I & II, we will then use the cluster sample to 
explore not only the orbital characteristics of the individual 
clusters (Paper III), but global properties of the Galactic disk 
(Paper IV), including: (1) the local Galactic rotation curve 
and velocity field near the Sun, (2) the kinematics of the disk 
across the frontier separating R <Ro and R > Rq, and (3) the 
validity of the assumption of Galactic dynamical symmetry 
(e.g., north vs. south, Galactic quadrants I/II vs. IV/III). 

In §6 and Table 12 of this paper we present the derived 3D 
space motions of the clusters that enable these future contribu- 
tions. In the preceding sections of this paper we explain how 
we selected our target clusters (§2.1 and §4) and which stars 
within each cluster field to probe (§2.2), the spectroscopic 
observations and the derivation of radial velocities (§3), and 
the means by which membership within each cluster is estab- 
hshed (§5). 

2. SOURCE SELECTION 

2.1. Cluster Sample Selection 

Our selection of specific open clusters starts with the 205 
clusters explored in the Dias et al. (2001, 2002a) catalogs, 

' The WIYN Observatory is ajoint facility of the University of Wisconsin- 
Madison, Indiana University, Yale University, and the National Optical As- 
tronomy Observatories. 
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which derive cluster membership using the statistical method 
of Sanders (1977). We also adopt the following criteria: (1) 
the clusters must have at least ten stars with Tycho-2 proper 
motions in the fields selected by Dias et al. (2001, 2002a), and 
(2) the cluster diameters cannot be much larger than the Hydra 
field of view (40'- CTIO, 60'- WIYN) so that the cluster can 
be sampled with a significant number of fibers. In addition, 
to obtain the greatest leverage on the local Galactic rotation 
curve the selected clusters span a wide area over the Galactic 
^gc-Ygc plane and reach to a heUocentric distance of >2.5 kpc. 
Neither age, distance from the Galactic plane, nor metallicity 
was considered as a selection criterion. 

Table 1 shows the basic cluster parameters of our sample 
with data taken from the Dias et al. (2002b) catalog, includ- 
ing coordinates of right ascension and declination (cols. 2 
and 3) and Galactic longitude and latitude (cols. 4 and 5), 
heliocentric distance (col. 6), log(age/years) and visual di- 
ameter of the cluster in arcminutes (cols. 7 and 8), and the 
observing run on which the cluster was observed (see be- 
low and Table 3 for definitions). The Galactic distribution 
of our final cluster sample of 71 clusters is shown in Figure 
3. The smaller number of clusters we have sampled in the 
I = 0-180° half of the Galaxy is result of a smaller amount 
of observing time obtained for the WIYN observations; how- 
ever future work in the research program wiU aim to remedy 
this deficiency. Figure 3 also shows the distribution of the se- 
lected cluster ages and distances from the Galactic midplane 
as a function of their Galactic radius (assuming the Sun is at 
8.5 kpc). More than half of our final sample have ages less 
than 200 Myr but older than 10 Myr (Table 1). The large 
number of relatively young clusters is important for kinemati- 
cal studies of the Galactic disk because, in general, open clus- 
ters should develop increasing deviations from "normal" disk 
rotation due to the scattering by molecular clouds over time 
(Spitzer & Schwarzschild 1951, 1953); however, clusters that 
are too young may still reflect the specific dynamical envi- 
ronment of their birth and may not yet have circular orbits 
(Lynga & Palous 1987). Figure 3c shows that all clusters in 
our sample are within 500 pc of the Galactic plane, and most 
are within 200 pc. This supports the notion that most clusters 
in our sample are likely to be "well behaved" in the sense that 
they have not been scattered far from the Galactic midplane 
and therefore are likely to still be on near-circular orbits (of 
course, we will revisit this question when we examine cluster 
orbits in detail, in a future contribution). 

Less than 25% of our clusters have estimated metalUci- 
ties ([Fe/H]), so we have little leverage on this aspect of our 
sample; however, we hope to derive metallicity estimates for 
some of our clusters in the future, using not only improved 
isochrone fits to CMDs aided with our membership data, but 
the spectra themselves. 

2.2. Stellar Selection Within Each Cluster 

Given that bulk 3D motions are our primary goal, the first 
stellar targets within each cluster selected for observation 
were those Hipparcos and Tycho-2 stars used in the Dias et 
al. (2001, 2002a) survey for the clusters. Because constraints 
on fiber optic placement with the Hydra instrument (i.e., two 
fibers cannot be closer than 25" in the Hydra setup) mean that 
in some cases not all desired stars can be observed in a cluster 
field, we must prioritize stars within a fiber setup. For this 
reason, stars were ranked in priority order based on the Dias 
et al. (2001, 2002a) derived membership probabilities, from 
highest to lowest probabiHty. Dias et al. (2001, 2002a) de- 



rived these probabilities based on the proper motions using 
the method of Sanders (1977). 

Next, additional Tycho-2 stars available in the Hydra field 
of view, but not used in the Dias et al. (2001, 2002a) study 
(because they lie beyond the cluster radius studied by these 
authors) were added as the next priority to the target list. 
For the WIYN/Hydra runs, no targets beyond the Tycho- 
2 stars needed to be selected because the combination of 
a smaller number of available Hydra fibers (90 vs. 132 
for CTIO/Hydra) and larger field of view (60' vs. 40' for 
CTIO/Hydra) typically meant that nearly all target fibers were 
filled with Tycho-2 stars. 

For the CTIO runs and for fields having less than 50 stars 
with available Tycho-2 proper motion data, we selected at 
lowest priority two additional sets of stars; first, stars be- 
tween V = 13-15 magnitude from the USNO-Bl.O catalog 
from within the cluster radius (with that value taken from the 
Dias et al. 2002b catalog), with the goal of searching for ad- 
ditional cluster members fainter than the V ~ 13.5 magnitude 
limit of the Tycho-2 survey, and second, we aUowed unused 
"field orientation probe stars" (FOPS; USNO Bl.O stars with 
12 < /?2 < 13) to be added to the bottom of the target priority 
list. At either WIYN or CTIO, fibers that were not assigned 
to targets were used for sky observations, with at least six 
(WIYN) or ten (CTIO) fibers positioned on random sky for 
sky subtraction of the stellar spectra. 

3. SPECTROSCOPIC SURVEY DATA 

3.1. Spectroscopic Observations 

We have collected homogeneous spectroscopic observa- 
tions for 71 open clusters using the HYDRA multi-fiber 
spectrographs on the Blanco 4-m telescope at Cerro Tololo 
Inter- American Observatory (CTIO) and the 3.5-m WIYN^ 
telescope at Kitt Peak National Observatory (KPNO). This 
project was conducted using publicly competed NOAO time 
and was granted long-term status^, which permitted observa- 
tions in the semesters 2002A-2004A over a total of fourteen 
awarded CTIO 4-m and six WIYN 3.5-m nights. The data for 
the clusters listed in Table 1 were obtained the nights of UT 

2002 March 8-12 ("Run 1"), 2003 March 16-21 ("Run 2"), 

2003 July 20-23 ("Run 3"), and 2003 August 2-8 ("Run 4") 
from CTIO. For more efficient observing, some cluster obser- 
vations scheduled for two August 2003 CTIO nights were in- 
terspersed with other targets for two other observing projects 
awarded telescope time over the course of eleven CTIO/Hydra 
nights in July and August 2003 (Runs 3 and 4). The WIYN 
data were observed on the nights of 2003 September 14-18 
("Run 5"). 

The CTIO observations made use of 132 Hydra fibers 
that are simultaneously dispersed onto a 2048 x 4096 pixel, 
SITe400mm CCD using the 380 grating with 1200 lines mm"' 
and with the fiber ends viewed by the spectrograph through 
the 100/um sht plate to improve the resolution to a dispersion 
of 0.68 A per resolution element (R ~ 15,000). The spectral 
range covered was 7740-8740 A. Data obtained at WIYN dis- 
persed the 90 Hydra fibers dispersed onto a 2048 x 2048 pixel 
CCD in the Red Bench Camera using the echelle (316@63.4) 
grating in 6th order; this yielded a dispersion of 0.82 A per 
resolution element {R ^ 13, 000) which was centered on the 

* The WIYN Observatory is a joint facility of the University of Wisconsin- 
Madison, Indiana University, Yale University, and the National Optical As- 
tronomy Observatories. 

^ This project was selected as an NOAO Ph.D. thesis project for PMF. 
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8220-8800 A spectral range. Typical signal-to-noise ratios 
of 10 or better were obtained; all cluster stars presented 
have at least S/N > 5. To aid the RV cahbration, multiple 
RV standards were observed on each run, where each "ob- 
servation" of an RV standard entails sending the Ught of the 
calibrator down 2-12 different fibers, yielding many dozen in- 
dividual spectra of each RV standard. We present here the 
results from analysis of spectra for 3436 individual stars out 
of 3537 with sufficient S/N observed in the fields of 71 open 
clusters. ^ 

3.2. Data Reduction 

Preliminary processing of the two-dimensional data was 
undertaken using standard IRAF^ techniques as described 
in the IRAF ccdproc documentation. After completing 
the CCD bias subtraction, overscan correction and trimming, 
the two-dimensional images were corrected for pixel-to-pixel 
sensitivity variations and chip cosmetics by applying "milky 
flats" according to the prescription outUned in the CTIO Hy- 
dra manual by N. Suntzeff.'" 

After basic processing the data were run through the IRAF 
routine dohydra. One dimensional spectra for each star 
were extracted from the two-dimensional CCD images and 
wavelength caUbrated with respect to a comparison lamp 
spectrum. Exposures of the PENRAY (CTIO; He, Ne, Ar, 
and Xe) or CuAr (WlYN) lamps were taken at each Hy- 
dra pointing through all fibers to provide comparison spectra 
yielding at least 1 1 prominent emission Unes roughly evenly 
distributed over the observed wavelength range. These com- 
parison spectra provide a wavelength solution (i.e., pixel to 
wavelength conversion) for each extracted object spectrum. 

3.3. Stellar Radial Velocities: Standard Stars 

All radial velocities were derived using IRAF's fxcor 
package, which we used first to determine RVs for the stan- 
dard stars. Radial velocity standard stars are used to check for 
systematics in the data, to determine the measured RV preci- 
sion, and to calibrate the zero-point of the velocity scale. The 
reduction to radial velocities employed essentially the classi- 
cal cross-correlation methodology of Tonry & Davis (1979). 
The template star input to the correlator is prepared from a 
high signal-to-noise (S/N) standard star spectrum from which 
the stellar continuum is fitted and subtracted. The result- 
ing spectrum is high and low pass Fourier-filtered to remove 
both high frequency noise (e.g., cosmic rays) and the low fre- 
quency variation cause by difference in instrument through- 
put. 

We first measured the RVs of standard stars by cross- 
correlating each Fourier-filtered standard star spectrum 
against every other standard star spectrum from its corre- 
sponding observing run over the restricted wavelength range 
of 8220-8680 A, which avoided possible contamination from 
nearby atmospheric fines. The resulting RVs from different 
cross-correlations for each individual standard star spectrum 
were averaged and the standard deviations measured; the re- 
sults are presented in Table 2. The average velocity stan- 

Of the 3537 stars observed with 5/^ > 5, 101 were peculiar stars (e.g., 
carbon stars. Be stars, young emission-line stars) that are excluded from the 
RV analysis using "normal star" cross-correlation templates discussed here. 

' IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of IJniversities for Research in As- 
tronomy, Inc., under cooperative agreement with the National Science Foun- 
dation. 

http://www.ctio.noao.edu/spectrographs/hydra/hydra-nickmanual.html 



dard deviation for the individual Hydra standard star spectra 
is (Ty < 2 km s"' for spectra with S/N > 20. 

We determined the level of the random and any unknown 
systematic RV errors from a prescription described in Vogt 
et al. (1995), which is based on the analysis of repeatedly 
observed stars (Table 2). In this case these stars were typ- 
ically observed through different Hydra fibers. The Tonry- 
Davis Ratio (TDR; Tonry & Davis 1979) for each spectrum is 
measured using fxcor. Since the TDR scales approximately 
with S/N we can, foUowing the method described in Vogt et 
al. (1995), determine approximate Icr errors in the RVs corre- 
sponding to a given TDR as: 



error Vr ■ 



(1-1- TDR) 



(1) 



The parameter a is a constant calibrated by the standard star 
data using the following formula, which is predicated on the 
assumption that the TDR is a good measure of the relative 
S/N, and where autocorrelations are not included: 



E,E/l+TDR/.;)-(V^..,-(V^7))^ 

X50,n 



(2) 



where Vr.ij is ith observation of the jth standard star, and (Vr.j) 
is the mean RV of the Jth standard star. We obtained the val- 
ues of Xso for our sample's number of degrees of freedom, 
where x^o n critical value of the distribution at the 

50% confidence level multiplied by n degrees of freedom, as 
described fully in Vogt et al. (1995). 

Since our target stars were selected based on proper motion 
criteria, they span a wide range of spectral types, including 
anything from hot O and Be stars to cool carbon stars. As 
a result, we observed a range of RV standard star templates. 
However, due to the lack of International Astronomical Union 
(lAU) RV standards hotter than spectral type AO, we used B 
and A stars from Fekel (1999) to provide RV standards for our 
hot star spectra. For a better match in the spectral types be- 
tween targets and cross-correlation templates, we divided the 
observed standard stars into "red" or "blue" subsamples for 
our cross-correlation templates. Stars were considered "red" 
stars if the Ca n infrared triplet (8498A, 8542A, and 8662A) 
was present in the spectra; this encompasses cool F through 
early M type stars. "Blue" stars have dominant Paschen se- 
ries lines and virtuaUy no Ca II triplet (i.e., O through A type 
stars). 

Initially we anticipated using primarily late type stars for 
our analysis, so that in the earlier runs (March 2002, March 
2003 and July 2003 observations; Table 1 - Runs 1, 2, 3) 
we did not obtain "blue" standards. Later it became clear 
that good RVs for hotter stars could be derived and we be- 
gan to collect blue standards. To cover the lack of blue stan- 
dards in the earlier runs, the "blue" August 2003 standards 
were used to reduce all CTIO "blue" target stars. This ap- 
proach was adopted and found to work moderately well be- 
cause even across observing runs all spectra were taken with 
the same instrument setup, are dispersion corrected uniformly, 
and should experience no flexure problems because Hydra 
uses a bench-mounted spectrograph. In the end we did find 
some offsets in the RV zero-points for some of the runs — 
but, ironically these were for the runs where we actually did 
take blue standards (see Sections 3.6 and 3.7). 

The standard star spectra (Table 2) provide a data set for 
calibration of the RV errors for each run according to Equa- 
tions 1 and 2. Table 2 lists each observation of a standard star. 
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the UT date of the observation (col. 2), its spectral type (col. 
3), which observation frame and fiber were used (cols. 4 and 
5), the TDR (col. 6), the mean derived radial velocity (Vr) and 
the measured standard deviation (cols. 7 and 8), followed by 
the average Vr and standard deviation for the cluster (cols. 9 
and 10), also the lAU or Fekel (1999) RV is shown for com- 
parison in col. 11. Using the Vogt et al. (1995) technique and 
measurements from Table 2, we find the parameters given in 
Table 3, which are used to determine errors for each of the 
runs. The variation in the a values can be due to the effects of 
focus and small spectrograph setup variations. 

3.4. Radial Velocity Standard Verification 

To check the reliability of our measurements of radial ve- 
locity standard stars, we compare the difference between our 
measurements (Table 4) and the lAU or Fekel (1999) values, 
as shown in Figure 4. As may be seen, the derived RVs for 
the standard stars are all within 2.3 km s~^ of the LAU values. 
We find that the difference between our measurement and the 
lAU values are no more than 2 times larger than the quadra- 
ture errors (as shown in Table 4); however, we find that the 
differences are randomly distributed and that the mean offset 
is less than 1 km s~' . Therefore we find there are no system- 
atic trends between our measured velocities and the cataloged 
values for the lAU standards, though we did find an offset for 
the "blue" Fekel (1999) stars, a situation that is analyzed in 
more depth below (§3.6). 

3.5. Stellar Radial Velocities: Target Stars 

Target stars were analyzed using the same reduction pro- 
cedure as the standard stars, with the exception that stars 
that showed both the Paschen series lines and any hint of 
the Ca II triplet were considered "green" stars. As with the 
RV standards, the targets were sorted into "red," "blue," and 
"green" sub-samples based on visual inspection of their spec- 
tral features in order to match them to the appropriate cross- 
correlation template. Figure 5a shows the 2MASS color 
distribution of the stars selected for each sub-sample. The 
"green" stars were tested against both templates to find the 
best match. Nearly all "green" stars became part of the "red" 
sample. 

Each group of target stars was processed through IRAF's 
f xcor package to cross-correlate them against the standard 
of the corresponding color class for their respective observing 
run, as described in §3.3. "Green" stars were cross-correlated 
against both red and blue templates and the derived RV was 
taken from the template that provided the better result. The 
final 2MASS color distribution for stars fitted with the "red" 
and "blue" templates is shown in Figure 5b. Uncertainties for 
the stars fitted to the "red" or "blue" templates were deter- 
mined using the a values from Table 3 and Equation 1 . 

3.6. Internal Comparison (Red vs. Blue) 

As an additional check on the measured RVs, we tested 
the internal consistency delivered by the separate "red" and 
"blue" reductions for a given run. To do this, stars in our 
"green" sample were correlated with both the "red" and 
"blue" standards. For stars with measured uncertainties in 
their blue measurement of less than 10 km s"' , Figure 6 shows 
iVr,biue-Vr,red) VS. Vr.red, where Stars with the best blue uncer- 
tainties (< 6 km s"') are denoted with black squares. The 
subsample of stars with uncertainties of < 6 km s"' was then 
fitted with a Une to determine if any zero point offset was 



needed between the red and blue samples. The fit to the data 
in the Vr^red vs. Vr^biue plane is given by: 

Vr^adopted = flQ + «! * Vr^blue- (3) 

The resulting fits for the March 2002, 2003, August 2003 and 
September 2003 runs are given in Table 5. We did not find a 
systematic difference between the < 10 km s"' and <6 km s"^ 
samples, just a larger scatter for the < 10 km s"' data. 

We find that the blue data are offset from the red data by a 
significant amount for the August 2003 and September 2003 
runs, and a small offset is found for the March 2003 run. To 
verify that it is the red values that are more reliable and to 
support the rationale that we shift the blue system RVs to the 
red system, we next compare our RV measurements from the 
red sample with those of previously pubUshed, high RV reso- 
lution studies for a number of open clusters. 

3.7. Systematic Effects and Comparison to Previous Results 

As an additional test of the reliability of our RVs, we have 
found previously published values for stars in nine clusters 
that we have observed. We undertook a comparison of our 
velocity measures to those in the following studies: IC 4561 
(Mermilliod et al. 1995; Meibom, Andersen, & Nordstrom 
2002), IC 4756 (Mermilliod & Mayor 1990), NGC 2099 
(MermiUiod et al. 1996), NGC 2423 (MermiUiod & Mayor 
1990), NGC 2447 and NGC 2539 (MermiUiod & Mayor 
1989), NGC 2682 (M67; Mathieu et al. 1986), NGC 5822 
(Mermilliod & Mayor 1990), and NGC 6134 (Claria & Mer- 
milliod 1992). In Table 6, we present a direct star-by-star 
comparison of RV results to the seminal work on M67 by 
Mathieu et al. (1986). Table 7 provides star-by-star com- 
parison for the other clusters listed above, which include the 
Tycho-2 star name, the star name from the corresponding pho- 
tometry reference used for identification in the previous RV 
studies, the stellar coordinates, our RV (Vr) and its uncer- 
tainty, the reference RV and its uncertainty, and the per star 
difference in these measurements. For stars in common be- 
tween the surveys, we find overall excellent agreement in the 
determined per star kinematics. 

In Figure 7, we compare differences between our own red 
data and previously published RVs as a function of photomet- 
ric parameters (e.g., magnitude and color) of the stars, where 
the data are color-coded by observing run (red = March 2002, 
green = March 2003, cyan = August 2003, and blue = Septem- 
ber 2003). We find no systematic trend with magnitude or 
color as shown in Figure 7a-d. The bottom two plots in Fig- 
ure 7 (panels e and f) show a comparison of AVr vs. Vr- While 
there may seem to be an odd trend at -25 > V, > -40 km s"' 
in panel (e), this is due mainly to one cluster — IC 465 1 — 
that has a peculiar offset. This is demonstrated by the "dis- 
appearance" of that odd trend when IC 465 1 is removed from 
the distribution (Figure 7f ; see §6.2.2 for a detailed discussion 
oflC 4651). 

Therefore, we find that our "red" sample, which we have 
made our standard reference, is consistent with previous work 
and this bears out our having corrected the blue sample RVs to 
the red RV system. The cause of this offset it probably due to 
a combination of using "blue" standards from different runs, 
as well as the fact that the two blue Fekel standard stars have 
only a few good lines for RV determination combined with 
large rotations with both Fekel stars having Vsin? 18 km 
s"'. 

4. FINAL CLUSTER SAMPLE 
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Table 8 summarizes all clusters observed, including UT 
date of the observation and exposure times, the numbers of 
stars selected to be cluster members by Dias et al. (2001, 
2002a) that were targeted with Hydra fibers (col. 4), the total 
number of stars and number of Tycho-2 stars observed (col. 

5) , the number of total observed stars with reliable RVs (col. 

6) , and the Tycho-2 stars with reliable RVs (col. 7). For the 
WIYN data, we were able to observe nearly 75-80% Tycho-2 
stars used in the corresponding Dias et al. (2001, 2002a) sur- 
vey. For the CTIO runs, we found that we were generally able 
to observe 50-80% of the Dias et al. selected Tycho-2 stars, 
and, in addition, sample an average of ^ 50 more non-Tycho- 
2 stars (since the latter were generally fainter by 1-3 magni- 
tudes, a lower fraction of them delivered reliable velocities in 
the allotted observing time). Since we are obtaining data for 
most of the Dias et al. stars, we will be able to compare our 
new membership data directly against the membership analy- 
sis done by these authors (see §6.2.3). 

5. CLUSTER MEMBERSHIP ANALYSIS 

One of the most complicated problems affecting studies of 
open clusters is membership contamination associated with 
their location within the densely populated Galactic plane. 
Large numbers of disk stars unrelated to the cluster lie along 
the CMD sequences of the typical open cluster and, given 
the typical motions of many objects within the Galactic 
plane, usually with rather similar velocities. To determine 
the bulk motion of clusters one must first isolate true clus- 
ter members from the dominant field star population in the 
fore/background. To accomplish this discrimination we have 
modified a previously implemented method designed to do 
just that. The star's proper motion, RV, and spatial distribu- 
tion are all used as inputs for a kernel-based, probabiUty dis- 
tribution function technique, described below, that eventually 
allows the cluster bulk motion to be determined from stars 
with high membership probabilities. 

5.1. Non-Parametric Frequency Function 

To determine cluster membership probabilities for stars 
based on RV and proper motion, we have chosen to use an 
empirical, non-parametric technique — modified from that 
described in Galadi-Enriquez, Jordi, & TruUols (1998) — 
that incorporates a kernel estimator (Hand 1982) to isolate the 
phase space distribution of cluster stars in a field. 

While we adopt the basic technique used by Galadi- 
Enriquez et al. (1998) for proper motions alone, we have gen- 
eralized it also to operate on a spatially-constrained, 1-D RV 
distribution as well as an RV-constrained, 2-D proper motion 
distribution. In principle, one could use either distribution 
separately for culling cluster members, but for the most se- 
cure assessment of membership we depend on the joint prob- 
abiUty distributions. This means, therefore, that we can only 
use stars having both RV and proper motion data. To improve 
our results further, we remove stars with large measurement 
errors in either proper motion or RV, or those stars that clearly 
have halo-like RVs, using the following constraints applied to 
the data: 



• II error limit: ^/cr 



+ cr^ < 10 mas yr 



The modified version of the Galadi-Enriquez et al. (1998) for- 
mulation is intended to perform better for our particular sur- 
vey circumstances — i.e., fewer numbers of stars per cluster, 
but high-quality RV data for these stars. Throughout the fol- 
lowing description we will demonstrate the basic features of 
our analysis via the example processing of the cluster NGC 
2682 (M67), for which the raw data are shown in Figure 8. 

5.2. 1-D Kinematical Distribution: Radial Velocities 

For our data, the RV distribution is found to be the most 
sensitive discriminator of cluster membership because of the 
small measured relative RV errors. When applying a ker- 
nel density estimator the empirical density function (V"^/) is 
comprised of both the cluster (c) and the field (/), where here 
V stands for the RV distribution. Since the observed empiri- 
cal density function is the sum of two underlying distributions 
(e.g., ipl^f = one must decompose the distributions 

to isolate the cluster function. Because of the accuracy of the 
RV data and the small intrinsic velocity dispersions of open 
clusters (0.5-3 km s"'), we expect to be able to discriminate 
the cluster and field fairly readily. To do so, however we must 
first isolate the field population to verify which peak in the 
ip^^j distribution is due to the cluster Differences in the clus- 
ter versus field distribution should be evident by looking at 
samples of stars drawn from different radii from the cluster 
center. A useful initial assumption is that stars outside of the 
cluster radius are "non-members", and these can provide a 
reasonable estimate of ipj. 

The RV data kernel analysis is comprised of four steps: (1) 
All RV data are convolved with a Gaussian kernel to homoge- 
nize our errors for a given cluster. This kernel has a width de- 
termined by the mean RV errors from all of the observed stars 
in a given cluster field. Because open clusters have intrinsic 
velocity dispersions of 1-3 km s~' in addition to our measure- 
ment errors, we Umit the Gaussian width to be at least 3 km 
s"' and at most 10 km s~'. Applying the Gaussian kernel to 
smooth our RV data (fpc+f^ produces the smoothed field plus 
cluster distribution an example for NGC 2682 (M67) is 
shown in Figure 9a. (2) We apply the same Gaussian kernel 
to smooth the RV data of stars that are outside the cluster ra- 
dius (utilizing the cluster diameters from Dias et al. 2002b). 
This smoothed RV distribution is used as the field distribution 
^'^(Figure 9b). (3) We wish to determine the probability of 
any particular star with a given RV being a member of the 
cluster, so we need to determine the normalized probability 
distribution: 



^^^/V„)-^j^(V„) 



(4) 



RV Umit: -200 < K < 200 km s"' 
RV error Umit: cry^ < 10 km s~^ 



The cluster probabiUty distribution is shown Figure 9c); 
however, we see that a few outUers, which are non-member 
stars within the cluster radius, are still visible in the distri- 
bution. (4) We assume that the strongest peak in the "clus- 
ter" probability distribution, P^ , belongs to the cluster, and 
perform a 1-D Gaussian fit to this peak (Figure 9d; dotted 
line). This Gaussian fit is used to determine RV membership 
probabilities, Pj , for all stars in the field and to exclude non- 
member RVs that still may appear in '^^ . 

5.3. 2-D Kinematical Plane: Proper Motions 

The proper motion kernel analysis is similarly comprised 
of four steps, but now applied in 2-D. This technique for 
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proper motions is identical to that used in Galadi-Enriquez 
et al. (1998) with the exception that instead of using a spatial 
membership separation (which we used for the RV distribu- 
tion described above) to establish we have chosen to use 
the RV separation described above (i.e., those stars outside 
the Gaussian fit to the RV distribution, Pj = 0, are considered 
the "field" population). The proper motion kernel uses the 
following equation analogous to that used in the RV analysis 
above: 



(5) 



where a' is acos(S). Continuing our example of NGC 2682, 
we apply the 2-D kernel smoothing to the proper motion dis- 
tribution as shown in Figure lOa-d. The Gaussian fit to the 
field-subtracted distribution is used to determine proper mo- 
tion membership probabilities (Figure lOd) for stars in 
each cluster Both the RV and proper motion kernel analy- 
sis was performed on all clusters in Table 6. 

5.4. Calibration of the Membership Criteria 

To determine the RV membership "cutoff" criteria, we have 
chosen to analyze in detail one of our best sampled clus- 
ters, our example NGC 2682. Using techniques standard 
dynamical techiques from Pryor & Meylan (1993), we per- 
formed an iterative 3a rejection using the full sample of NGC 
2682 RV data. We find an intrinsic velocity dispersion of 
(Tin, = 0.96 ± 0.29 km s"' . As a comparison, using proper mo- 
tions, Girard et al. (1989) found that NGC 2682 (M67) has 
aint = 0.81 ± 0.10 km s~'. Comparing the RV member stars 
left after applying the iterative 3a rejection, we find that all of 
the remaining stars have Pj > 70%. More lenience is given 
to the proper motions due to the larger average error, and in 
this regard we follow the criterion used by Dias et al. (2001, 
2002a). As a result, we have chosen to define cluster mem- 
bership as stars that have Pj > 70% and > 5 1%. 

5.5. Results of the Membership Analysis 

Cluster membership was determined by jointly assessing 
the probabilities from the 1-D RV distribution (Pj) and the 2- 
D proper motion distribution (P^). The probabilities for each 
star analyzed in the cluster NGC 2682 are included along with 
the RV and proper motion data in Table 9. This table includes 
the star name from the Tycho-2 survey, or if not a Tycho- 
2 star, another identifier (for M67 we have IDs from Eggen 
& Sandage 1964; Sanders 1977; Montgomery, Marschall, & 
Janes 1993; Fan et al. 1996). The table then lists, in order, 
the right ascension and declination for each M67 star (cols. 
2 and 3), the Tycho-2 proper motions and errors (cols. 4-7), 
our measured RV and error (cols. 8 and 9), and which spec- 
tral cross-correlation template was used to derive these (col. 
10). In addition, we have included the membership probabil- 
ity from Dias et al. (2001, 2002a; col. 11) for comparison 
to our derived membership probabiUties P^ x 100 (PM; col. 
12), X 100 (RV; col. 13), and P'J", the joint probabiUty 
(P'f = pVpK X 100; col. 14). The stars selected as cluster 
members are presented in boldface type. 

Similar probability data are given for the other clusters in 
our sample in Table 10, which is available in electronic for- 
mat. In this table we give for each star observed its Tycho- 
2 name, or, if a non-Tycho star, an identifier with the format 
"XXXX_f_####" for the added "filler" candidate FOPS guide 
stars or "XXXX_u_####", for USNO B-1.0 catalog "filler" 



stars, categories described in the observational criteria in §2.2. 
In §6.2.3 our analysis of the cluster memberships of these 
stars are compared against the membership analysis by Dias 
et al. (2001), whose membership probabilities are based only 
on proper motion. 

5.6. Cluster Membership and Cluster CMDs 

As shown in Figure 1 , with only photometric data the iden- 
tification of open cluster sequences in the CMD can often be a 
tricky prospect. Our radial velocity cluster memberships can 
significantly aid in clarifying the location of these cluster se- 
quences. The 2MASS and Tycho-2 photometry for all stars in 
our survey with measured RVs are listed in Table 1 1 . 

Figure 1 1 shows the 2MASS CMD for the example cluster 
NGC 2682 with our spectroscopically-observed stars identi- 
fied, and with large circles denoting stars selected to be mem- 
bers based on both RV and proper motion. Triangles denote 
stars that have P^ > 70% but which do not have Tycho-2 
proper motion data. For now we present CMDs without red- 
dening corrections applied, because this is a non-trivial pro- 
cess in that not all hne-of-sight reddening (the values typically 
given in catalogs such as Schlegel, Finkbeiner, & Davis 1998) 
is necessarily foreground to the cluster. One can see from the 
CMD that in this case our membership census yields mem- 
bers that fall primarily along the photometric sequences of 
M67 apparent in the CMD. Similar 2MASS CMD member- 
ship plots for all clusters we have studied are shown in Fig- 
ures 12-15. As in the case of M67, our identified members 
typically fall in the expected locations of the main sequence 
turn-off (MSTO) or giant branches of the clusters, when those 
are obvious; however, in many cases the CMDs are crowded 
with field star contamination and our identified members help 
clarify the cluster sequences. This is particularly useful in the 
fairly common situation where the giant branches are sparsely 
populated. As we shall show in another contribution (Frinch- 
aboy et al., in prep), our ability to clarify the CMD locations 
of cluster giant branches and MSTOs greatly improves the 
isochrone fitting for these systems. 

6. KINEMATICAL RESULTS 

6.1. Derived Cluster Space Velocities 

The cluster bulk RV is calculated using cluster members 
(e.g., as shown in Table 9) and techniques from Pryor & Mey- 
lan (1993) to determine the cluster mean RV and error in the 
mean. The cluster mean bulk proper motions are calculated 
using the following equations (and a symmetrical version for 



^ Macos((5) 



En 
i=l 



(6) 



(7) 



The derived cluster bulk motions are given in Table 12, 
where we fist the numbers of members with full space mo- 
tions (col. 2) and the 3D members plus the stars determined to 
be members by RV criteria alone (3D+RV; col. 3), along with 
the resulting bulk kinematics and the associated uncertainties 
(RV from all 3D members; col. 4), RV from 3D and additional 
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"RV only" members (col. 5), and equatorial and Galactic sys- 
tem proper motions (cols. 6-9). We find two clusters NGC 
1513 and NGC 7654 with only one star selected for member- 
ship (i.e., the membership method found no more than one 
star with a given RV within the errors); given the uncertainty 
in selecting among single star subsamples to define the actual 
"cluster", we remove these two clusters from further analysis. 

6.2. Comparison to Previous Results 
6.2.1. NGC2682 (M67) Example 

In § 3 .7 and Table 6 we have already demonstrated a star-by- 
star comparison of derived RVs for the example cluster M67. 
For stars in common between the surveys, we find excellent 
agreement in the determined per star RVs (previously shown 
in Figure 7). Now we compare the derived bulk space velocity 
for this very well-studied cluster to the most detailed, previous 
studies of M67. 

In Table 13 we compare our derived mean proper motion 
and radial velocity for M67, averaged over these measured 
piffameters for 10 stars we determined to be reliable 3D mem- 
bers of the cluster, against derivations of these bulk motion 
parameters by other authors. With regard to to the previously 
derived bulk RV for M67, our mean radial velocity is consis- 
tent with previous measurements by Mathieu et al. (1986) and 
Scott et al. (1995), and lies within 0.2 km s~' of the rather pre- 
cise value given in the Mathieu et al. study. The total number 
of published clusters having as extensive and detailed RV cov- 
erage as the Mathieu et al. M67 study is less than ten, whereas 
our study now provides high precision RVs for stars in nearly 
five times as many clusters. We also find proper motion re- 
sults more or less consistent with previous measured values, 
with our ^5 value being bracketed by the measurements 
by Dias et al. (2001) and Kharchenko et al. (2005) results and 
our cos 5 reasonably close to the values for this proper mo- 
tion component derived by these two other studies. The previ- 
ous studies have smaller errors in their mean due to the larger 
numbers of "member" stars used in the determination of the 
bulk proper motion. 

Thus we find that our survey results are consistent with the 
very detailed analysis of previous M67 work. Despite the fact 
that M67 is probably one of the most well-studied open clus- 
ters in the Galaxy and previous studies typically utiUzed many 
more stars than we have, our results deliver comparable pre- 
cision to the best of these because of the greater purity of our 
samples, and, in the case of the RV measurement, the velocity 
resolution of our spectra. 

6.2.2. Comparison to Previously Derived Bulk Cluster 
Radial Velocities 

A compilation of our derived mean cluster RVs compared 
to those found previously by other authors is given in Table 
14. We have found previous results for 25 of our 71 studied 
clusters, some with multiple studies. In general, we find con- 
sistency with the previous studies to the few km s~' level as 
shown in Figure 16. but in a few cases, there are more sub- 
stantial differences. 

Figure 16 shows that the clusters NGC 457, NGC 884, and 
NGC 957 have discrepant RVs found between our work and 
any previous study; however all of these clusters, plus NGC 
2264, were studied by Liu, Janes, & Bania (1989). The Liu 
et al. (1989) study is comprised of only a few possible cluster 
members observed (e.g., for NGC 884 and NGC 957 only two 
stars each and these clusters also have large mean errors). We 



believe our results, which incorporate both RV and proper mo- 
tion membership, are superior to those from Liu et al. (1989). 
Even with the small numbers of stars in both studies, we find 
that our results are marginally consistent with Liu et al. (1989) 
for NGC 2264. 

6.2.3. Comparison to Previously Derived Bulk Cluster 
Proper Motions 

In Table 15, we compare our derived open cluster bulk 

proper motions with the previous results of Dias et al. (2001, 
2002a). The latter surveys used only the Tycho-2 proper mo- 
tions to derive membership and the cluster bulk proper mo- 
tions. Table 15 compares the numbers of stars used by Dias et 
al. and their derived mean cluster proper motions (col. 2-5) 
to our own sample statistics and derived mean proper motions 
(col. 6-9). As shown in Figure 17 (grey histogram), three 
clusters — Collinder 258, Lynga 1, and NGC 6250 — show 
large inconsistencies (A^ > 5 mas yr"') between our results 
and those of Dias et al. We also reminder the reader that we 
have already excluded two other cases (NGC 1513 and NGC 
7654; see §6.1) from our study, because we identified only 
one star selected as a possible cluster member. Looking fur- 
ther at the proper motion difference outliers, we find that each 
Lynga 1 and NGC 6250 have only one star with fully derived 
3D kinematics and in the case of ColUnder 258 there are only 
two member stars. Thus, we conclude that our analysis may 
have settled on the wrong star(s) to represent the cluster in 
these cases and that the results for Collinder 258, Lynga 1, 
and NGC 6250 (in addition to NGC 1513 and NGC 7654) 
may not be reliable. For the remaining 66 of our 71 clusters, 
our "re-measured" proper motions are within the Icr errors of 
those found by Dias et al. (2001, 2002a), though our data gen- 
erally have comparable or smaller resulting errors in the mean 
(as shown in Figure 17) of ~ 1.5 mas yr~^ 

The direct comparison to the Dias et al. proper motions is 
shown in Figure 17, with the full sample shown in grey and 
various subsamples based on the number of members in either 
survey shown by the colored histograms. A somewhat close 
agreement with Dias et al. is expected because we are deriv- 
ing proper motions using a subsample of Dias et al. stars and 
adopting the same astrometry. A key difference, however, is 
that a number of Dias et al. "member" stars are excluded by 
our RV membership criterion so that, while we typically de- 
rive approximately the same bulk motions as Dias et al. these 
authors allow many more actual non-members to enter their 
sample; nevertheless, that Dias et al. include more actual non- 
members seems to have relatively small effects because these 
authors are typically averaging over a large number of stars 
in each cluster, including, apparently, sufficient numbers of 
true members to get close to the correct proper motion. We 
show in Table 15 the numbers of Dias et al. member stars 
{Poias > 50%) that are confirmed to be members (col. 12) and 
how many we find to be unlikely members (col. 13) based 
on the addition of our RV analysis. On average we find half 
of the Dias et al. "member" stars to be non-members when 
we account for the RVs. This suggests that use of proper mo- 
tion data of the quality of Tycho-2 alone may be insufficient 
to determine reliable cluster memberships, though, when av- 
eraged over many multiple stars and applying the 3o- rejec- 
tion of outlier proper motions adopted by Dias et al., these 
proper motions are useful for deriving the cluster bulk proper 
motion. The Dias et al. membership inaccuracies are likely 
lessened for closer clusters (e.g., d <2 kpc) which have more 
bright Tycho-2 stars. Using the sub-samples from Figure 17, 
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we see that when both samples have a lot of "members" there 
is convergence to a common proper motion, as expected. We 
also see that as the sample sizes decrease the measured proper 
motion differences grow. It is clear that, at least in our case, 
when we have too few stars we may have trouble "finding" 
the true cluster members (e.g., as in the examples of NGC 
1513 and NGC 7654). However as both our and the Dias et 
al. also studies drop to a few stars per clusters, it is difficult 
to determine which study is correct. We argue that given our 
more restrictive 3D membership criteria that ours is superior, 
though further study will be needed to confirm this assertion. 
Thus, while Tycho-2 has the best currently available astro- 
metric data, more strict RV discrimination such as we provide 
can substantially improve the application of these data for de- 
termining cluster motions, given a sufficient number of RV 
members. 

7. SUMMARY 

We have derived high precision (typically < 3 km s~' un- 
certainties) radial velocities for 3436 stars in the fields of 71 
open clusters within 3 kpc of the Sun. This represents the 
largest sample of clusters assembled thus far having uniformly 
determined, high-precision radial velocities. To extend this 
uniformity to the other velocity dimensions, our survey has 
focused primarily on obtaining spectra of stars having mea- 
sured Tycho-2 proper motions; however, our target list was 
appended with other stars in the cluster fields to expand the 
membership census for each cluster We have jointly applied 
three criteria — spatial position, radial velocity and proper 
motion (in two dimensions) — to derive high quaUty clus- 
ter membership probabilities for the samples stars. In at least 
half of our clusters we have found at least three stars in the 
field that are reliable members of the cluster using all of these 
criteria. 

Using these member lists, we have averaged the RVs and 
the Tycho-2 proper motions to derive mean space velocities 
for each cluster With few exceptions, our mean cluster RVs 
are close to those previously derived for the several dozen 
clusters that have been surveyed by other groups. A compari- 
son of our mean cluster proper motions with those by Dias et 
al. (2001, 2002a) — who also relied on Tycho-2 proper mo- 
tions — shows that both data sets are in general agreement, 
though our results should be more rehable given our more 
stringent assessment of cluster membership (i.e., we add high 
quaUty RVs to the proper motion criteria used by Dias et al.). 
We find that typically a large fraction of the Dias et al. stars in 
each cluster field do not meet our most restrictive, joint mem- 
bership criteria. In a few cases with discrepant proper motion 
results compared to those derived by Dias et al. we find that 
the differences may be due to a critically small numbers of 
stars surviving our 3D "membership" criteria; i.e. in some 
of these cases (namely Collinder 258, Lynga 1, NGC 1513, 
NGC 6250, and NGC 7654) it is likely that our results, based 
on only one or two stars, might be wrong due to the improper 



identification of cluster members. Nevertheless, our data pro- 
vide reliable 3-D space motions for 66 open clusters. 

In most cluster fields we have explored, our membership 
analysis provides valuable new benchmarks for improved 
isochrone fitting of the cluster CMDs, which is useful for 
estimating ages, distances, metallicities and/or reddenings to 
these systems. The resulting distances and metallicities will 
allow a new attempt at measuring the Galactic metallicity gra- 
dient with these clusters. With improved distances and more 
reliable space velocities, the orbits of the clusters can be de- 
rived under an assumed Galactic potential and solar Galacto- 
centric distance. Alternatively, these space velocities can be 
used as tracers of the local velocity field and be used to inves- 
tigate the Galactic rotation curve with a set of objects having 
velocity independent distances and uniformly derived, quality 
space velocities. We intend to address these science issues in 
future contributions in this series. 

Finally, our census of reliable cluster members provides a 
primary target list for future efforts to explore these open clus- 
ters with either high resolution spectroscopy or high precision 
astrometry, like that expected from SIM PlanetQuest. 
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0, 


.2126 


1033 


O A1 A 

/.olO 


7 


2 


NGC 2670 


08: 


■ A C 

.4j, 


:30 


/I Q 

— 4o, 


■ AQ 

.4o, 


:00 


262.1476 


0, 


.7868 


1 188 


O AQA 

/.oyu 


7 


2 


NGC 2669 


AQ. 

Uo, 


■ /I A. 
.40, 


:22 


CO. 

— j2. 


■ CA 

.JO, 


:54 


267.4854 


—J, 


AOCA 

.o2jU 


1 A/1 A 

1U4d 


O OOO 

/.y2 / 


OA 

2U 


2 


Trumpler 10 


08: 


■ An. 

:4/, 


:54 


—42: 


:27: 


:00 


262.7906 


0, 


.6740 


424 


7.542 


29 


1 


NGC 2682 


AC 
\JO. 


Ol, 


:18 


+11, 


.45, 


:00 


122.9232 


oo 
—LI . 


A/1 AQ 

.U4UU 


QQQ 

yu5 


Q /1QQ 

y.4uy 


oc 
2D 


1 


Collinder 205 


AOi 

uy, 


■An, 
■UU, 


:32 


AQ 

— 4o, 


• CO. 

.jy 


:00 


269.2091 


1 

— 1, 


QAIA 


1 QCO 

15Dj 


O OQA 

/.2UU 


c 
D 


2 


IC 2488 


AOi 


■2/, 


:38 


CO 

— D / 


.AA 

.UU 


:00 


277.8298 


A 

—4, 


A 1 QO 

.4iy2 


1 1 0/1 

1134 


5.1 13 


1 

15 


1 


NGC 2925 


AO. 


■ ■3 "2. 

.■53: 


:11 


CO 

— JJ 


.oo 


:54 


274.6855 


1, 


■7570 


774 


O OCA 

/.5D0 


1 A 

lu 


1 


NGC 3680 


11; 


■ OC 

;2j; 


:38 


A1 

— 4j, 


.1/1 
:14, 


:36 


124.9390 


— 1, 


OOOA 

.LLLKi 


QOO 

935 


9.077 


D 


2 


Collinder 258 


12: 


■ OO. 

;2/: 


:10 


A A 

— OU 


./I A 

:4o. 


:00 


299.9710 


1, 


QAC/1 

,yOD4 


1 1 0/1 

1 lo4 


O AOO 

O.032 


D 


1,2 


NGC 5281 


1 1. 

I J. 


.40, 


:35 


AO 
—02, 


■ CC. 

.J J, 


:00 


309.0102 


—2, 


/1Q1 C 

.4y 1 J 


1 1 AQ 
1 lUo 


O 1 /I A 

/. 14d 


/ 


2 


NGC 5316 


13: 


:53: 


:57 


— 61 : 


:52: 


:06 


311.6017 


2, 


. 1 144 


1215 


8.202 


1 /I 

14 


2 


Lynga 1 


14, 


■ AA. 

■ UO, 


:02 


AO 

— o2. 


:09: 


:00 


310.8493 


— 0, 


.3373 


2283 


8.007 


3 


2 


NGC 5460 


14, 


■ AO. 


:27 


AQ 

— 4o, 


:20: 


:36 


316.3148 


5, 


.6067 


678 


8.207 


35 


2 


Lynga 2 


14, 


:24: 


:35 


A 1 
— Ol , 


:20: 


:00 


313.8642 


— 0, 


.4544 


1000 


8. 122 


10 


1 


NGC 5617 


1 A 
14, 


■ OA 

.2y , 


:44 


AA 
— OU, 


■ /10 
.42, 


:42 


317.5264 


2, 


ACC 1 

.UoDl 


1 coo 
1D33 


O Q1 c 

/.yiD 


1 A 
lU 


2 


NGC 5662 


14: 


:35: 


:37 


—56: 


:37: 


:06 


319.5288 


4, 


.5444 


AAA 

oOo 


O QAO 

/.yoo 


29 


1 


NGC 5822 


1 ^ 

1 J, 


■ A/l ■ 
.U4, 


:21 


C/l ■ 
— j4. 


■ OO. 

.2j, 


:48 


324.3610 


1 , 


oo A 1 

. / 2U 1 


Q1 O 

yi / 


Q QO 1 

o.o21 


oc 
3D 


2 


NGC 5823 


15: 


:()5: 


:30 


—55: 


:36: 


:12 


343.8165 


19, 


.8092 


1 1 AO 

1 iy2 


QAA 

o.you 


1 o 
12 


2 


NGC 6025 


1 A. 
lo. 


■ AQ. 

■Uj. 


:17 


A A 

— OU 


■ OC. 

.2j, 


:54 


329.7454 


—2, 


OA/1 Q 

.2U4o 


OCA 

/JO 


O QQQ 

/.ooy 


1 /I 
14 


1 


NGC 6031 


lo. 


■ AO. 

:U/. 


:35 


C/1 , 

— j4. 


■ AA 

:0U 


:54 


327.7257 


— D, 


/10CA 

,42D0 


1823 


O AAA 

o.ooy 


3 


2 


NGC 6067 


Id; 


■ 1 ■i. 
.15: 


:11 


C/1 

— j4. 


. 1 o 


:06 


127.7404 


2, 


AOOQ 

.Uo/U 


1 /I 1 o 

141 / 


O AOA 
5.U /O 


1 A 

14 


2 


Harvard 10 


10, 


■ 1 Q, 
.10, 


:48 


C/1 

— j4. 


.CA. 

.DO 


:00 


329.8356 


■3 

—J, 


OQ/1/1 
.2044 


1 1 o 

1312 


0/1A 

5.34U 


oc 
2D 


2 


NGC 6124 


10, 


■o<, 
.2D, 


:20 


ylA 

— 4U 


.jy 


:12 


332.9179 


■7 
—J, 


1 AAQ 

.1000 


CIO 

D12 


Q 1 /IO 

5.14/ 


OQ 

3y 


1 


NGC 6134 


10; 


■OO, 

:2/; 


:46 


ACt 

—49 


.no 

:uy 


:06 


335.2223 


1 

— 1, 


/1000 

.42/2 


Ql O 

913 


O QAO 

5.905 





2 


Ruprecht 119 


10, 


.OQ, 
.2(5, 


:15 


CI . 

— jl 


.OA 

.jU 


:00 


333.2758 


1 

— 1, 


QOQ/1 

.o /y4 


QCA 

yDO 


A QCO 

0.5D3 




5 


1 


NGC 6167 


10; 


;j4; 


;34 


ACt 

—49 


:4d. 


:18 


338.4047 


1 

1, 


O 1 AA 

,2100 


1 1 QO 

llUo 


7.887 


7 


1 


NGC 6250 


1 A' 

10, 


■ <o. 
.J /. 


:56 


/I c 

— 4j, 


■ CA 

.JO, 


:12 


341.9974 


1 

—1 , 


C 1 AA 

.D loO 


OAC 

oOD 


O /I 1 c 
/.41D 


1 A 


4 


NGC 6281 


i / , 


.U4. 


:41 


—J / , 


■CO 


:06 


345.2791 


—J, 


AC A/1 
.Uj04 


/10Q 

4 /y 


Q /1QO 

o.4y / 


Q 
O 


2 


IC 4651 


17: 


:24: 


:49 


-49: 


:56: 


:00 


340.0881 


-7, 


.9068 


888 


9.057 


10 


1 


NGC 6405 


17: 


:40: 


:20 


-32: 


:15: 


:12 


356.9316 


-1, 


.5491 


487 


7.974 


20 


2 


NGC 6416 


17: 


:44: 


:19 


-32: 


:21 


:42 


357.9402 


-1 


.6054 


741 


8.087 


14 


2 


NCG 6603 


18: 


:18: 


:26 


-18: 


:24: 


:24 


15.8996 


0, 


.3505 


3600 


8.300 


6 


2,3 


IC 4756 


18: 


:39: 


:00 


+05: 


:27: 


:00 


36.3807 


5, 


.2422 


484 


8.699 


39 


4 


NGC 6705 


18: 


:51: 


:05 


-06: 


:16: 


:12 


15.3951 


-9, 


.5927 


1877 


8.302 


13 


4 


NGC 6811 


19: 


:37: 


:17 


+46: 


:23: 


:18 


73.9778 


8, 


.4808 


1215 


8.799 


14 


5 


NGC 6866 


20: 


:03: 


:55 


+44: 


:09: 


:30 


60.3897 


-6, 


.0501 


1450 


8.576 


14 


5 


NGC 6885 


20: 


:11: 


:58 


+26: 


:29: 


:00 


66.1352 


-6, 


.3113 


597 


9.160 


10 


4 


Berkeley 86 


20: 


:20; 


:24 


+38: 


:42: 


:00 


76.6667 


1, 


,2725 


1112 


7.116 


6 


5 


Platais 1 


21: 


:30; 


:02 


+48: 


:58: 


:36 


92.5613 


-1 


■6461 


1268 


8.244 


10 


5 


NGC 7209 


22: 


:05; 


:07 


+46: 


:29: 


:00 


102.7010 


0, 


.7820 


1168 


8.617 


14 


5 


NGC 7654 


23: 


;24: 


:48 


+61 


:35: 


:36 


117.2878 


10, 


.8044 


1421 


7.764 


15 


5 



"Run 1: March 2002 (CTIO), Run 2: March 2003 (CTIO), Run 3: June 2003 (CTIO), Run 4: August 2003 (CTIO), Run 5: 
September 2003 (WIYN) 



Open Clusters as Galactic Disk Tracers I 



TABLE 2 
RV Standard Observations 







Avg. 




lAU 


UTDate Star Spec. Type Frame* Fiber TDR (kms-i) 










Full table given in journal 



^Due to poor S/N this star was not used as a cross-corelation template 
"Not an lAU RV standard star 



TABLE 3 
Calibration of RV Errors 







Star 


Degrees of 








Rim Dates 


Ruii# 


Colors 


Freedom 




<\ 


Teleseoiie 


2002 Mar 10-13 


1 


Red 


24 


23.337 


35.19 


CTIO 4-m 


2003 Mar 16-21 


2 


Red 


22 


21.337 


27.69 


CTIO 4-m 


2003 July 19-22 


3 


Red 


41 


40.334 


50.24 


CTIO 4-m 


2003 Aug 01-07 


4 


Red 


72 


71.333 


33.59 


CTIO 4-m 


2003 Aug 01-07 


4 


Blue'' 


9 


8.343 


84.87 


CTIO 4-m 


2003 Sep 14-17 


5 


Red 


6 


5.3481 


53.22 


WIYN 3.5-m 


2003 Sep 14-17 


5 


Blue 


6 


5.3481 


41.14 


WIYN 3.5-m 



"These data were used as "blue" cross-correlation templates for all CTIO runs. 



TABLE 4 

Comparison of RV Standard Star Values 



Star 


Spectral Type 


UT date 


V, (km s"' ) 


(km s ' ) 


Differenee (km s ') 


HD 


126053 


GOV 


11 Mar 2002 


-20.8 ± 0.6 


-18.5 ± 0.4 


-2.3 ± 0.7 


HD 


150798 


K2n-in 


11 Mar 2002 


-1.8 ±0.8 


-3.7 ± 0.2 


+1.9 ±0.8 


HD 


157457 


G8III 


12 Mar 2002 


+17.7 ± 0.6 


+17.4 ±0.2 


+0.3 ± 0.7 


HD 


136202 


F8III-IV 


22 Mar 2003 


+54.9 ± 0.9 


+53.5 ± 0.2 


+1.4 ±0.9 


HD 


157457 


G8III 


19 Mar 2003 


+17.5 ± 0.5 


+17.4 ± 0.2 


+0.1 ±0.5 


HD 


168454 


K2.5IIIa 


19 Mar 2003 


-20.7 ± 0.4 


-20.0 ± 0.0 


-0.7 ± 0.4 


HD 


9138 


K4III 


21 Jul 2003 


-34.2 ± 0.6 


-35.4 ± 0.5 


+1.2 ±0.8 


HD 


18884 


M1.51IIa 


21 Jul 2003 


-25.7 ± 1.5 


-25.8 ±0.1 


+0.1 ± 1.5 


HD 


9138 


K4111 


22 Jul 2003 


-34.2 ± 0.3 


-35.4 ± 0.5 


+1.2 ±0.6 


HD 


18884 


M 1.5111a 


22 Jul 2003 


-25.8 ± 0.7 


-25.8 ±0.1 


-0.0 ± 0.7 


HD 


107328 


K0.5IIIb 


23 Jul 2003 


+37.4 ± 0.5 


+35.7 ± 0.3 


-1.7 ±0.6 


HD 


146051 


M0.5III 


23 Jul 2003 


-18.8 ± 0.8 


-19.8 ± 0.0 


+1.0 ±0.8 


HD 


693 


F5V 


02 Aug 2003 


+14.1 ± 1.4 


+14.7 ± 0.2 


-0.6 ± 1.4 


HD 


9138 


K4III 


02 Aug 2003 


-34.6 ± 0.6 


-35.4 ± 0.5 


+0.8 ±0.8 


HD 


18884 


M 1.5111a 


02 Aug 2003 


-25.7 ± 0.7 


-25.8 ±0.1 


+0.1 ±0.7 


HD 


693 


F5V 


03 Aug 2003 


+14.8 ± 1.6 


+14.7 ± 0.2 


+0.1 ± 1.6 


HD 


18884 


M1.5IIIa 


03 Aug 2003 


-25.5 ± 0.5 


-25.8 ±0.1 


+0.3 ± 0.5 


HD 


693 


F5V 


05 Aug 2003 


+14.4 ± 1.0 


+14.7 ± 0.2 


-0.3 ± 1.0 


HD 


693 


F5V 


06 Aug 2003 


+ 14.8 ± 2.1 


+14.7 ± 0.2 


+0.1 ± 2.1 


HR 


7773 


B9IV 


06 Aug 2003 


+1.2 ±2.8 


-1.0 ±0.2" 


+2.2 ± 2.8 


HR 


675 


A2V 


07 Aug 2003 


-1.6 ±4.5 


+0.4 ± 0.2" 


-2.0 ± 4.5 


HD 


18884 


M1.5IIIa 


07 Aug 2003 


-24.6 ± 1.4 


-25.8 ±0.1 


+1.2 ± 1.4 


HD 


693 


F5V 


08 Aug 2003 


+14.8 ± 1.2 


+14.7 ± 0.2 


+0.1 ± 1.2 


HD 


9138 


K4III 


08 Aug 2003 


-34.7 ± 1.2 


-35.4 ± 0.5 


+0.7 ± 1.3 


HD 


18884 


M1.5IIIa 


08 Aug 2003 


-25.8 ± 1.2 


-25.8 ±0.1 


+0.0 ± 1.2 



"Stars HR 675 and HR 7773 are not lAU standards. Due to the lack of lAU standards hotter than F type stars, we 
used stars from Fekel (1999). 



TABLE 5 

Comparison OF Measured Red vs. Blue RVs 



Run 


ao 


ai 


RMS 


Blue Correction 










(lans-i) 


March 2002 


-0.47 


0.97 


3.06 


0.0 


March 2003 


-1.26 


0.99 


3.74 


-1.3 


August 2003 


-10.09 


1.00 


4.43 


-10.0 


September 2003 


-6.05 


0.98 


4.63 


-6.0 



Frinchaboy & Majewski 



TABLE 6 

NGC 2682 Comparison TO Previous RV Results 













"2000 




^2000 


Vr 




VrMathwii 




AVr 


Star (Tycho-2) 


SAND 


FBC 


MMJ 


ES 


(hi) 




(°) 






Ikms-') 


,kms-') 


(kms-1) 


08 1 ^-01 1 

i ^ \J 1 i ^ 1 1 


0721 


2554 


6475 




8-SfV4Q (SS 


I [ 


35 '08 9 


-t-33 01 


79 


4-34 2 


9 


—1 2 




0727 


2425 


5158 




8-5n-44 98 


1 1 


37:30.4 


4-2 1 .44 


1-82 


4-16-3 




5.1 


yjO It- VJi. J J 11 


0978 


3364 


6482 




8-'=; 1 -1 7 4.8 


I [ 


45'22 7 


-t-35 17 


59 


4-34 7 


0-6 


5 


\JO It — \J U J 11 


1010 


3536 


6485 


1008 


8-'^ 1 ■?? 81 

O.J 1 O 1 


\ [ 


48"01 8 


-t-33 00 


0-69 


4-33-6 


0-4 


-0.6 




1016 


^'^4Q 


6484 




O.J i . i / . 11^ 


Y y 


48' 16 2 


-t-35 66 


v.jy 




n 7 

U- / 


1 4 


OSM-OlfS^l-l 
in — \7 1 U J i 1 


1023 


3668 


6487 


1001 


O.J 1 .Z.VJ. Ot 


I [ 


48-40 5 


-t-3 54 


1 01 


4-3-6 


1 


—0 1 


0814.-01 1 lQ-1 

\JO It — \^ 1 i 17 1 


1040 


3566 


6488 




O.J 1 .i. J . / o 


I [ 


49-49 4 


-t-34 60 


82 


4-32 


6-3 


2.6 " 


0814-01763-1 


1054 


3347 


6489 


4020 


8:51:17.03 


11 


50:46.4 


4-32.63 


0-72 


4-33-5 


0-4 


-0.9 


0814-01099-1 


1074 


3204 


6492 




8:51:12.70 


11 


52:42.4 


4-34.50 


0-68 


4-34-1 


3-1 


0.4 


0814-01823-1 


1221 


4149 


6497 




8:51:43.56 


11 


44:26.4 


4-36.73 


0-58 


4-34-1 


3-1 


2.6" 


0814-01515-1 


1250 


3755 


6499 




8:51:29.91 


11 


47:16.8 


4-31-06 


0-99 


4-34.3 


1.4 


-3.2 


0814-01007-1 


1277 


4117 


6502 


3032 


8:51:42.32 


11 


50:07.8 


4-33.18 


0.78 


4-34.0 


0.5 


-0.8 


0814-01147-1 


1279 


3726 


6503 




8:51:28.99 


11 


50:33.1 


4-33.93 


0.66 


4-33.3 


0.4 


0.6 


0814-01471-1 


1288 


4118 


6505 


3034 


8:51:42.36 


11 


51:23.1 


4-32.25 


0.67 


4-33.3 


0.4 


-1.0 


0814-01225-1 


1293 


4039 


6050 


3035 


8:51:39.38 


11 


51:45.4 


4-33.86 


0.97 


4-34.1 


0.5 


-0.2 


0814-00795-1 


1306 






3065 


8:51:49.36 


11 


53:38.9 


-2.44 


0.63 


-1.1 


0.4 


-1.3 


0814-00134-1 


1402 


4878 


6508 


2152 


8:52:10.97 


11 


31:49.2 


4-33.12 


0.76 


4-33.6 


0.4 


-0.5 


0814-00847-1 


1327 


3979 


6507 


3086 


8:51:37.18 


11 


59:02.4 


4-11.62 


0.87 


4-12.0 


0.4 


-0.4 


0814-02313-1 


1585 


5191 






8:52:26.33 


11 


41:27.7 


4-33.25 


1.31 


4-34.1 


0.1 


-0.9 



References. — SAND: Sanders (1977), FBC: Fan et al. (1996). MMJ: Montgomery et al. (1993), ES: Eggen & Sandage (1964), Mathieu: Mathieu et al. (1986). 



Star listed at spectroscopic binary in Mathieu et al. (1986). 



TABLE 7 

Star-b y-Star Comparison to Previous Results 



0:2000 
Star (hr) 


52000 
(°) 


Vr Vr (Other) AVr Other Other 
(kms"') (kms"') (kms"') ID Ref. 


Full table given in journal 


References. — 1 


Mermilliod et 


a. (1995), 2: Meibom et al. (2002), 3: Mermilliod & Mayor (1990), 4: Mermilliod 



et al. (1996), 5: Mermilliod & Mayor (1989), 6: Claria & Mermilliod (1992) 
^Star listed at spectroscopic binary. 



Open Clusters as Galactic Disk Tracers I 



TABLE 8 

Statistics of Observed Cluster Stars 



Cluster 


UT-Date 


Exposures 




# Observed" 


# Vr Stars 










(total/TYC) 


(total/TYC) 


NGC 129 


900^ Spn 1 5 


1 V 600 e 


29 


67/ 


29 


29/29 


NOP '^Rl 

iNw^-. JOl 


700^ Spn 1 5 


1 V 600 V 


20 


65/ 


20 


20/20 


NOP d'>l 


700^ Spn 17 


1 V 600 e 


17 


61/ 


17 


17/17 


i>VJ\^ OOt 


700^ Spn 1 S 


1 V 600 e 


23 


45/ 


23 


23/23 




700^ Spn 1 6 


1 V 600 ^■ 


19 


63/ 


19 


1 Q/l Q 




700^ Spn 1 5 


1 V 600 c 


19 


48/ 


19 


19/19 




700^ Spn 17 


1 V 600 V 

J A Ul.7\7i> 


43 


43/ 


43 


43/43 




700^ Ana 07 


1 V 600 e 

J A UUUil 


13 


108/ 


13 


70/13 


O HJUJv O 


700^ Spn 1 6 


1 V 600 e 

J A U1.7\7l> 


14 


64/ 


14 


14/14 




7001 "spn 1 7 


1 V 600 c 




37/ 


35 


'^5rt5 


NOP 7nQQ 


7001 Spn 1 S 


1 V 600 c 

J A UUUj 


37 


51/ 


37 


37/37 


If liiii*/^li**tilrr» 1 
iVlltUdldllvU 1 


zAfyJD OCp J. O 


1 y 600? 


37 


49/ 


37 


37/37 




900^ Mar 1 8 


■3 V 6OO7 


12 


55/ 


17 


22/13 


MfX' 9964 


9009 Mar 1 9 

ZfWZi IVXOI J.Zi 


1 y 600? 


18 


38/ 


34 


32/21 




900"^ Mar 1 8 


1 V 600c 

J A DUUA 




78/ 


00 


95/ 




9009 Mar 1 9 


■3 V- 600? 

D A UUUJ 


38 


95/ 


38 


57/38 


Nnr 9393 


9009 Mar 1 1 


1 V 600? 

D A UvUJ 


55 


96/ 


94 


63/40 


Nnr 9354 


7001 Mar 90 


1 V 600 c 


20 


107/ 


57 


78/52 


NPP 


7001 Mir 90 


1 V 600 c 


25 


100/ 


33 


58/20 


NPP 74.7'^ 


7009 Mqr 1 7 

ZI.A7Z IVltll IZ 


1 V QOOc 


50 


107/103 


96/63 


NPP 74.'^7 


7001 Mqr 1 Q 

Z.\}\)J IVltll 17 


1 V 600 c 

J A U 1.71.7 L> 


75 


107/ 90 


57/50 


NPP 74.47 


7001 Mnr 17 

ZI.A7J IVltll 1 / 


1 V 600 c 

J A Ul.7\7i> 


34 


81/ 50 


49/36 


NPP 74S9 


7009 Mir 1 1 

ZI.A7Z IVltll 1 J 


1 V 600 c 

J A U 1.71.7 1> 


32 


103/101 


67/29 


Nnr 2516 


2nn3 Mar 1 S 

Z-v/\y _J IVltll 1 O 


3 X 600 1 


45 


102/ 


37 


29/21 




7001 Mnr 90 
ZUU-J iVlttl ZU 


1 V 600 c 

D A OUUA 




117/ 





109/ 


Nnr 9597 


9001 Mar IQ 

ZUVJ J IVltll 1 7 


1 V 600 c 

J A U 1.71.7 1> 


30 


92/ 


48 


63/36 


NOP 7*^4*^ 


7009 Mar 1 1 

ZVJVZ IVltll 1 I 


1 V 600 c 

J A U 1.71.7 1) 


80 


102/ 


96 


54/49 


NOP 7'^47 


7001 Mar IQ 
zxjyjj ividi 1 y 


1 V 600 c 

J A UUUj 


19 


115/ 


30 


44/19 




0(\(W Mar 1 8 


1 V 600? 


'^n 
j\j 


111 


40 


41 /96 


T»jnP 9'%4.S 


9009 Mar 1 1 


a V QOO? 

D A 7UUi 




82/ 


81 


5Q/S4 




900^ Mar 1 Q 


1 y 600? 
D A UUUJ 


17 


81/ 


24 


55/16 


NnP 9'%70 


900^^ Mar 9 1 


1 y 600? 

J A DUUi 




lAI 


13 


'^6/1 n 

JO/ lU 


T»jnP 9^70 


900'^ Mar 90 


a V 600? 

J A OUUi 


Q 


56/ 


16 


96/ 7 


NOP 9fifi0 


900'^ Mar 1 8 

ZUU J IVlol 1 o 


1 y 600? 

D A UvUJ 


16 


100/ 


30 


64/21 


Trump ler 10 


7009 Mar 14 
ZUUZ IVldJ It 


1 y QOOp 

J A 7UUA 


22 


87/ 


32 


49/20 


NPP IfiRl 


7009 Mar 14 

ZI.A7Z IVltll 1*+ 


1 y 600 c 

J A UUUA 


28 


109/ 


28 


65/28 


Pnllinrlpr 905 


7001 Mar IQ 
ZI.A7 J IVltll 1 y 


1 y 600 c 


12 


67/ 


14 


32/ 9 


IC 2488 


7009 Mar 1 1 
Z1.7UZ IVlttl 1 I 


1 y 600 c 

J A UUUA 


40 


94/ 


93 


68/'^^ 

uo/o J 


NPP 9Q9'S 


7009 Mqr 1 7 

ZI.A7Z IVltll IZ 


1 y QOOc 


32 


98/ 


95 


71/44 


NOP '^(^sn 


7001 Mqr 17 

ZVA7 J IVltll 1 / 


1 y 600 c 


14 


63/ 


23 


47/23 


Pnllinrlpr 9'^S 


7009 Mqr 1 1 

ZI.A7Z IVltll 1 I 


1 y 600 c 


13 


92/ 


89 


77/74 




7001 Mqr 91 

ZW\7 J IVltll Z I 


1 y 600 c 




89/ 


20 


51/1 1 


KTPP S9S 1 


7001 Mir 1 Q 
ZUU J IVlttl 1 y 


1 y 600 c 


12 


81/ 


27 




NPP S^l 


7001 Mar 17 

ZI.A7J IVltll 1 / 


1 y 600 c 


25 


103/ 


77 


52/46 




7001 Mar IQ 

ZUU J IVldl 1 " 


1 y QOOc 

J A y\j\)!i 


9 


79/ 


22 


42/17 




9001 Mar 1 8 


1 y 600? 


40 


105/ 


63 


50/43 


T \mcj^ 9 
Liyiiga. z. 


9009 Mar 1 1 

ZVUZ IVlttl 1 J 


1 y QOO? 


13 


&!/ 


79 


62/15 


Nnp '>fii7 


9001 Mar 1 Q 

Z,UUO IVlol 17 


1 y 600? 

D A UUUj 


35 


88/ 


41 


42/23 


Nnp '>fifi9 


9009 Mar 14 

ZiUUZ IVXOl 1*+ 


■3 y QflO? 


60 


93/ 


84 


52/46 


Nnp '^899 


900'^ Mar 1 7 


1 y 600? 

D A UUUj 


140 


119/111 


70/68 


Nnp '^89'^ 


9001 Mar 1 8 

Z,UU,J IVlol 1 o 


1 y 600? 

D A UvUJ 


10 


100/ 


29 


57/23 


Nnp fi09S 


9009 Mar 1 1 
z\juz IVlol 1 1 


1 y 600? 

J A UUUJ 


66 


103/ 


92 


77/37 


Nnp ^ifi'^ 1 


9001 Mar 9 1 
ZUUO iVloI Zl 


1 y 600? 

J A UUUA 


\\ 


90/ 


21 


40/17 


Nnp 6067 


7001 Mar 1 7 
ividi 1 / 


1 y 600? 

J A UUUj 


24 


107/ 


77 


36/30 


± XCli- V CU- U 1 V/ 


2003 Mar 20 

^\J\J ^! IVltll 


1 V 600 c 


34 


118/ 


71 


49/40 


Nnr 61 94 

i > V.J \^ U 1 Z.'-r 


7009 Mnr 14 

Z1.7\7Z IVltll iT' 


3 X 900.V 


30 


91/ 


80 


59/53 


Nnr 6 1 34 


2003 Mar 1 8 

^ IVltll 1 O 


1 y 600 c 


60 


71/ 


23 


41/21 


NPP 6 1 «1 


7009 Mnr 1 4 

Zl.7l.7Z IVltll It 


3 X 900.V 


10 


86/ 


67 


48/ 8 


Riinrprht 1 1 Q 


7009 Mqr 1 1 

Z1.7V7Z IVltll 1 J 


3 X 9()0.v 


15 


96/ 


89 


77/22 


NGC 6250 


2003 Aug 05 


3 X 90().v 


14 


74/ 


21 


33/11 


NPP 69 S 1 


7001 Mqr 9 1 

Zl.7l.7J IVltll ZI 


1 y 600 c 

J A Ul.7l.7i) 


21 


81/ 


31 


47/22 


IC 465 1 


2009 Mar 1 3 

£^\>\}Z^ IVltll 1 J 


3 X 900.S- 


19 


96/ 


77 


78/30 


NGC 6405 


2003 Mar 21 


3 X 600.S- 


30 


110/ 


51 


70/50 


NGC 6416 


2003 Mar 21 


3 X 600i 


32 


100/ 


63 


52/34 


NGC 6603 


2003 Jul 21 


3 X 600i 


22 


110/ 


39 


55/23 


IC 4756 


2003 Aug 05 


3 X 600j 


56 


94/ 


74 


74/56 


NGC 6705 


2003 Aug 02 


3 X 600i 


32 


80/ 


46 


47/34 


NGC 6811 


2003 Sep 16 


3 X 600i 


64 


77/ 


77 


64/64 


NGC 6866 


2003 Sep 17 


3 X 600^ 


52 


62/ 


62 


52/52 


NGC 6885 


2003 Aug 04 


3 X 600i 


29 


110/ 


36 


57/29 


Berkeley 86 


2003 Sep 16 


3 X 600* 


46 


67/ 


67 


46/46 


Platais 1 


2003 Sep 18 


3 X 600^ 


48 


50/ 


00 


48/48 


NGC 7209 


2003 Sep 16 


3 X 600i 


42 


49/ 


49 


42/42 


NGC 7654 


2003 Sep 17 


3 X 600i 


14 


62/ 


19 


14/14 



^iVoio, is the number of Tycho-2 stars observed, which are taken fmm the study of Dias et al. 
(2001, 2002a). 
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TABLE 9 

Cluster Analysis - NGC 2682 Example' 





O!2000 




"2000 




fla COS 


fJ'6 


Vr 






Probabilities 




Star 


(hr) 


(°) 


(mas yr"') 


(mas yr""^) 


(kms"^) 


Type" 


Dias 


c 


pV 

c 


rtini 

n 


fQ7S4-MMTS 1 ■^9 

17 / J'T iVliVlJ J 1 JZ. 


O . Jl.* .T-Z.. JZj 


1 1:39:49.7 






J J . J J m J . Zt- 


Red 




0.0 


100.0 


0.0 




8-'S0-4'^ S7 

O . JVy.T' J.J / 




25 9 ± 1 7 


—32 7 ± 1 9 


—3 84 lb 1 07 


Red 





0.0 


0.0 





ns 1 ^ n9'^n9 i 


8-S0-44 08 

O. jU.'-f'+.70 


1 1 ■37'30 4 


—1 s -t- 9 9 


z.z ni Z.J 


9 1 44 -1- 1 89 

Zl .tt IC i .oZ 


Red 





1.3 


4.5 





n.o\}\)- r sj^^z.'^ I o 


8-S0-4fi 89 
o. jU. "-to. oZ 


1 1 ■52'57 4 






— 7(i 71 -1- S 47 
— / 0. / i zn j.t / 


Red 




0.0 


0.0 









XX m,JJ m\JOty 


—6.5 i 1.6 


-5.4 it 1.6 


^1 ni 4-0 79 


Red 


98 


87.9 


99.3 


87.3 


T4S99-FRr'9Rnn 


S-SO-SS 41 

O . J V../ . J O . 1 


1 1"58"14 4 






a 1 99 _L 4 an 

Jl.ZZ _l_ T^.Jl.7 


Red 




0.0 


88.7 







O.J 1 .IJ J. J 1 


11 "45 "02 7 


—8 8 it 1 9 


—7 2 it 1 8 


ao _|_ 9 04 


Red 


99 


72.1 


63.1 


45 5 




8-S1 -OS 89 
o. J i . Wj .oZ 


1 1 'SI' 10 7 






'in 1 -U 4 44 

J7.i7 ^ t.tt 


Red 




0.0 


53.3 





14794 FRr"^ni 7 


8 S1 -Ofi SO 
o. J i .LfO. JLf 


1 1 ■'^S-SO 






J / .DO ^ Z.yt 


Red 




0.0 


71.7 







8^1 -07 00 

O . J 1 .U / .7" 


I 1 ■'^8-9(S S 

I I .JO .ZU. J 






901 70 4- 9 6^ 
zwi. /y ^ z.uj 


Red 




0.0 


0.0 


0.0 


141 nQS FRT'^OQQ 


8-^1 -00 90 

O.J 1 .LI^.ZiLI 


1 1 ■^7-01 

1 1 . J / .\J 1 .K) 






37 16 zb 1 77 


Red 




0.0 


77.5 


0.0 




8-S1 ■ 1 9 9'^ 
o. J i . 1 z..iLj 


1 1 '47' 15 






J J. 7 1 zn 4-. J J 


Red 


c 


0.0 


89.6 





0S14-010QQ-1 




11:52:42.4 


—7.8 ± 2.0 


—4,5 it 2,0 


^4 cn -|_ /:8 


Red 


98 


90.0 


98.7 


88,8 




8S11'^ 10 

O. J 1 . U. IW 


1 1 ^701 
1 1 . J / .Ul .w 






■2a go _|_ oc 

JJ.7Z ^ Z.ZJ 


Red 




0.0 


99.8 


0.0 


UO It—Ul 7J 1-1 


8'Sl -14 '^S 

O. J 1 . It, J J 


1 1 ^S OO S 

1 1 .t J .UW. J 


—12.9 ± 1.8 


—12.8 it 1.8 


an ft4 _|_ 1 '28 

J7.Ut ^ I.JO 


Red 





0.0 


47.3 


0.0 


nsi4-ni448-i 


8:51:14.74 


1 1 '^O OQ Q 

1 1 . JW.U7.7 


36.5 ± 2.4 


ivJt.j ^ Z.J 


40 87 4- 76 

t7.0 / IE U. /U 


Red 




0.0 


0.0 


0.0 


UolM— Ul /DJ-1 


9.''>\ -17 0"^ 

O.J 1.1 / ,\JJ 


1 1 •'^0•4f^ 4 
1 1 . JW.tO.t 


—10 1 it 2 2 


— 6 S -1- 9 

— D.J ^ Z.J 


'29 f:'2 _|_ 79 

JZ.DJ ^ U. / Z 




QQ 

yy 


56.5 


97.9 


SS 4 
J J.t 


fiKid.nidO'^.i 


o<ox< X / >xu 


11 •48>1/i 9 


HI ^.u 


— rfS 8 -1- 1 Q 




ivcU 


OQ 


86.4 


92.0 


79,5 


VO J-"t J. ~ X 


0<*7X> X / t*tO 


11:45:22.7 


—7.3 it 1.8 


—5,7 it 1,7 


"^S 1 7 4- fl 


Red 


99 


98.5 


95.4 


94,0 


n^uu 1 ijv^jjv/-' 


8-S1 -91 8S 

0. J 1 .z, 1 .0 J 


1 1"43"17 8 






98 84 4- 3 08 

zo.ot- hi J. 70 


Red 




0.0 


61.9 


0.0 


141 OTQ-FRP^SOS 


8-^1 -91 0(^ 

O.J 1 .Z, 1 .7U 


1 1 -Si-oo 1 

1 1 .J J .V7 7. 1 






04 4- 4 40 

JT-.OT- m 'T.T^l./ 


Red 




0.0 


97.0 


0.0 


OKI 4-01 ^'^1-1 


Ji'«;i '99 81 


1 1 •48*ni 8 

X X . to. UX.o 


—7.9 it 2.0 


— ^ ^ -1- 1 Q 


00 -1- 


Red 


99 


98.9 


99.0 


07 Q 

y 1 my 


0S14-01 T 19-1 




11:49:49.4 


-6.6 it 2.5 


—4,7 it 2,9 


14 60 4- 89 


Red 


97 


83.9 


98.3 


82,5 


0814-01 01 1-1 


8-^1 ■9fi 4"^ 

O.J 1 .Z-U-T-J 


1 1 ■4'^-^0 7 

1 1 .T-J . JLI. / 


— 9 3 zt 1 8 


-6.2 it 1.8 


40 06 4- 9 17 
'-tKj.yxj IC z. 1 / 


Red 


99 


77.6 


33.0 


25.6 


0814-01^1^1-1 


8-^1 -If, 84 

O.J 1 .Z.U.Ot- 


1 1 ■48-40 S 

1 1 .T-O.tW. J 


— 16 4 zt 1 6 


— 6 1 zt 1 6 


3 54 zb 101 


Red 





0.0 


0.0 





\JO It — \J 1Z.VJ 1 


8-^1 -97 01 

0. J 1 .Z. / .W 1 


1 r5r57 6 

1 1 .J 1 . JZ.U 


—9 it 1 6 


—6.2 it 1.6 


49 93 zb 2 94 


Red 


99 


83.7 


0.0 





H94'?-FRr'571 fi 

1X7^*-' 1 ij V - ' / 1 vj 


8-51 -28 61 

O.J 1 .Z.O.VJ 1 


1 1 ■38-32 

1 l.JO.J^.L/ 






9^; A7 _|_ 1 nc 

Z-VJ.U/ 1 1.7J 


Red 




0.0 


36.0 


0.0 


141 91 S FRP'^790 


8-S1 -98 09 
0. J i .Lo.y L 


1 1 ■4'^-08 

i 1 .tJ.WO.7 






1 '1 -1- 4 44 
jZ.Ij ziz t.tt 


Red 




0.0 


95.4 



\}.\) 


0S14-01147-1 


K'Sl '98 


1 1 'SO'H 1 


—9.5 it 1.8 


—5,3 it 1.7 




Red 


99 


74.9 


99.8 


74,8 


0$i14-01^1^-1 


8:51:29,91 


11:47:16,8 




— /; 9 -1- 1 ^ 

D,^ ^ X,J 


11 -1- QQ 


IVcU 


yy 


94.8 


86.9 


82,4 


l41141-FRri780 
jn 1 jH- 1 -rDV_, J / 07 


8-Sl -10 80 

O.J 1 .JU.OW 


19 04 16 

lZ>.Wt. lU.W 






^9 9-2 _|_ 9 94. 

jz.zj m z.zt- 


Red 




0.0 


96.0 


0.0 


111 vJOL^O lUl 


8*S1*11 77 
0. J 1 . J 1 . / / 








19 88 4- 3 OS 

JZ.OO ^ J.UJ 


Red 




0.0 


98.7 


0.0 


1 1 1 ZttvJ XjiJ t.v/ 1 \J 


8'Sl 'X) 41 

O.J 1 . JZ-.T-l 


1 1 ^fi ^S 8 

1 1 .tw.tj.o 






'27 70 j_ 9 OQ 


Red 




0.0 


70.5 


0.0 


UO I't-U^UO / I 


8'Sl 'X) SQ 

O.J 1 ,jz,,jy 


11:48:52.1 


—11.5 it 1.8 


—5.0 it 1.8 


'2/; 07 4- 1 74 

JU.W / ^ l./t 


Red 


94 


27.7 


88.7 


24.5 


14190S FRP'^R'sfi 
nixuo rDV_,jo JD 


8*'^1 '^9 fiO 

O.J l.j.i.DW 


1 1 •49*0'i 

1 1 .t.i.UJ.W 






■27 ftA _|_ c Q7 


IXCU 


c 


0.0 


71.7 



w.w 


T4191fi-F*s9ni0 


8*S1 •'^'^ "^8 

O.J 1 .J J. JO 


11:45:59.9 






■2'2 Q£ _|_ 9 '2'2 
JJ.7U ^ Z.JJ 


Red 




0.0 


99.8 


0.0 


141 1 R4 FRP^0^7 


O.J 1 .J J. J J 


ir34'32 1 






07 _|_ A nq 

JJ.O/ ^ t.7J 


Red 


c 


0.0 


99.9 





0814-0084.7 1 


8-^1 -^7 1 8 
. J i . J / . i 


11"59"02 4 


_9q 9 + 90 

Z7.Z ZIZ Z.W 


—2.1 it 2.0 


11 69 -1- 87 
i 1 .\jz HI w.o / 


Red 





0.0 


0.0 





f07'^(^-FRri08S 


8"5r37 37 


1 9-01-1 8 

1 Z.V7 J .17.0 






9A Q*;, _|_ c 09 

Zt^.7U J_ J.JZ 


Red 




0.0 


20.7 





0814. 0I99S-1 


8-^1 -"^O '^8 

O.J i . J 7. JO 


1 l"51'4i 4 


—7 8 lb 2 2 


—9 1 it 2 3 


-3Q Qf. _i_ 07 
jj.ou _i_ \>.y / 


Red 


92 


34.5 


99.9 


34 5 


08 14 01 007 1 


S-S 1 -49 ^9 
O.J 1 .tZ. jZ 


1 1 -SfVfW 8 
1 1 . JVJ.V/ / .0 


—10 7 zt 2 3 


—4 2 zt 2 4 


1 « -!- 78 
J J. i zn u. / 


Red 


96 


38.6 


99.5 


38 4 


Ok1d-01 471-1 


8*^1 '49 
o.^x.t^. jn 


11:51:23.1 


—8.5 zb 2.2 


—4.7 zb 2.3 


-19 9 c _|_ ft /:7 
jz.zD HZ u.o / 


Red 


99 


88.2 


96.0 


84.6 


0S14-0189'^-1 

\}0 in — W 1 OZ. J 1 


8-S 1 -41 ^fS 

0. J 1 J . JU 


1 1 ■44-26 4 


— 6 2 ± 2 


— 8 5 zb 2 


7a _|_ S8 
JU. /J in 1.7. JO 


Red 


89 


39.6 


81.9 


32 4 


f07S0-FRP41 ^8 


S-S 1 -44 0(S 

O.J 1 .T'T^. 7U 


1 1 ■18^50 1 

1 I . JO .J 7. J 






-29 _|_ 1 AO 
jz.uj in i.t7 


Red 




0.0 


97.9 





0814-009(^4-1 


8-S 1 -dft 00 
O.J 1 .^yj.y/y 


1 1 ■1(S^ 1 8 8 
1 i . ju. 1 0.0 


-12.8 it 1.6 


-59.8 lb 1.6 


07 07 4- 75 
/ .u / ^ W. / J 


Red 





0.0 


0.0 





141 940 FS9049 


0. J i / .t J 


1 1 ■47^00 8 
i i .t / . W7.0 






49 QO -U 9 00 


Red 


c 


0.0 


17.2 





14191 1 FRP4901 


8-'^ 1 -48 48 
O.J i -to. to 


1 1 ■42'23 7 






'!<; 19 _|_ 9 47 
JJ. IZ ^ z.t / 


l\CU 


c 


0.0 


95.4 



w.w 


0814 09047 1 

Wo It- V/z.V/'t / i 


8-'^1 -48 f\d 

O.J i .tO.Dt 


1 1 ■40- 1 f, 

i l.t7. i J.U 


_i 1 /; -1- 1 

1 i.U ^ i.U 


U.J ^ i.U 


49 Q< _|_ 9 10 

tZ.OJ ^ Z.iW 


IxCU 


OS 

7J 


25.6 


17.2 


4 4 


T419S1-F'?9041 

1 1 1 Z. J 1 tjOZ,V/T- 1 


8-Sl -48 70 
O.J 1 .to. / w 


1 1 •47-'^S 8 

1 1 .t / . J J.O 






88 95 4- 9 74 

00. ZJ ^ z. /t 


Red 




0.0 


0.0 


0.0 


0814 0070^ 1 

WO It-WU / 7 J- i 


8'S1 -40 '^fy 
O.J i .t!7. JO 


1 1 •'%l-'^8 Q 
i 1. jj. J0.7 


_i 1 -2 _J_ Q 

— 1 i. J IjZ U.!7 


"^8-1-0 
J.O ^ u.y 


— 9 4J. 4- 61 

— Z.tt ^ U.DJ 


IVCU 


w 


0.0 


0.0 



w.w 


fQ77'^-S ANDI ^07 


8*^1 •'il 78 
0. J 1 . J 1 . / 


12:04:47.1 






■29 QQ _|_ 1 QfZ 


Red 




0.0 


99.0 


0.0 


141484 FRP44'^'^ 


8-^1 '^"^ "^0 

O.J i . J J. JW 


1 1 -^fii-M 
i 1. JO. i / .w 






11 08 4- 1 S7 
JI.70 ^ J.J / 


ivCU 


c 


0.0 


94.2 



w.w 


0814-099^*^-1 

wo 1*+ WZ.Z, J J 1 


8'*;i '^4 01 

O.J 1 . Jt,7 1 


1 1 •40-96 8 
1 i.tw.z.w.o 


—22.3 ±1.6 


-45.8 lb 1.6 


4Q ^0 4- 78 

t7. JW ^ W. / 


Red 





0.0 


0.0 


0.0 


0814-01011-1 


8:51:56.01 


11:51:26.6 


-7.2 =t 1.6 


-14.6 lb 1.5 


29.79 it 2.21 


Red 





0.0 


72.9 


0.0 


H1422-FBC4502 


8:51:56.10 


11:39:14.0 






27.05 it 5.68 


Red 


c 


0.0 


40.4 


0.0 


H1473-FBC4508 


8:51:56.24 


11:53:35.9 






35.59 it 2.92 


Red 


c 


0.0 


92.8 


0.0 


H1449-FBC4515 


8:51:56.61 


11:47:25.0 






33.59 it 3.47 


Red 


c 


0.0 


100.0 


0.0 


H1468-FBC4578 


8:51:58.65 


11:52:15.0 






38.84 it 2.83 


Red 


c 


0.0 


56.9 


0.0 


H1492-FBC4654 


8:52:01.59 


12:01:03.2 






17.73 ±4.11 


Red 


c 


0.0 


0.0 


0.0 


H1452-FBC4706 


8:52:03.51 


11:47:48.0 






34.78 it 3.23 


Red 


c 


0.0 


97.5 


0.0 


H1414-FBC4777 


8:52:06.37 


1 1:37:30.7 






0.42 it 2.26 


Red 


c 


0.0 


0.0 


0.0 


0814-00134-1 


8:52:10.97 


11:31:49.2 


-10.6 it 1.9 


-3.7 lb 1.9 


33.12 ±0.76 


Red 


94 


35.2 


99.5 


35.0 


H1426-FBC4887 


8:52:11.40 


11:40:32.0 






31.94 ± 3.44 


Red 


c 


0.0 


94.2 


0.0 


H1601-FBC5015 


8:52:16.90 


11:48:31.0 






29.86 ± 6.40 


Red 


c 


0.0 


74.0 


0.0 


0814-02313-1 


8:52:26.33 


11:41:27.7 


-13.7 lb 2.4 


-7.8 ± 2.6 


33.25 ± 1.31 


Red 


21 


2.9 


99.7 


2.9 



Boldfat^i: enlrii^s donolo slar sck'Llcd as ckislci iiiciiihcrs iisiii!; iIk' .>D ci ilci ia. 

Type of spectral template used in the cross-correlation lor RV determination, as described in §3.3. 

Star not used in EJias etal. (2001) membership analysis. 
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TABLE 10 

Cluster Analysis - All data^ 



02000 52000 tiaCosS fi^ Vr Probabilities 

Star (hr) (°) (masyr-') (masyr-') (kms"') Type'' Dias ^ 

Full table given in journal 



''Type of spectral templale used in Ihc cross-coirclalion for RV dclcrniinalion, as described in §3.3 
''Star not used in Dias etal. (2001) membership analysis. 
Star not analyzed due to large errors (see §5.1). 



TABLE 1 1 

2MASS AND Tycho-2 Photometry for Targeted Cluster Stars 



Star 2MASSID 02000 ^2000 Bt^Ao Vrycho 2MASS / 2MASS 2MASS ifj 3DMemb? 

Full table given in journal 



star not matehed to any 2MASS point source. 
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TABLE 12 
Derived Cluster Bulk Kinematics 





Members 


Bulk Vr3D 


Bulk Vr,3D+Rv 


f^a cos S 






iJ-b 


Cluster 


3D 3D+RV 


(k„,s-l) 


(kms->) 


(mas yr^^) 


(mas yr ' ) 


(mas yr"') 


„™,y,-l) 



Berkeley 86 


2 




-25.54 ±2.60 




-3.80 ± 1.27 


-4.58 ± 1.24 


-5.88 ± 1.25 


0.90 ± 1.26 


CoUinder 205 


1 


6 


27.40 ±4.07 


28.18± 1.52 


-4.00 ±2. 10 


8.00 ± 1.90 


-8.08 ± 1.90 


-3.84 ±2. 10 


CoUinder 258 


2 


6 


16.04 ±0.60 


14.90 ±0.93 


-0.90 ±2.62 


-5.00 ±2.40 


0.35 ±2.61 


-5.07 ± 2.41 


Harvard 10 


3 




-18.17± 1.94 




-3.88 ± 1.07 


-11.64± 1.11 


-11.99± 1.11 


2.58 ± 1.07 


IC 2488 


3 




-1.84± 1.52 




-7.73 ± 1.56 


8.93 ± 1.50 


-9.84 ± 1.50 


-6.55 ± 1.56 


IC 4651 


9 


10 


-33.34 ±0.67 


-33.30 ±0.67 


-1.72±0.70 


-2.76 ±0.69 


-2.54 ±0.69 


2.04 ±0.70 


IC 4756 


7 


13 


-24.44 ± 0.75 


-25.08 ±0.64 


0.68 ±0.74 


-1.99 ±0.72 


-1.95 ±0.72 


-0.76 ±0.74 


Ktiarchenko 1 


3 




12.74± 1.31 




6 10± 1 23 


-1.91 ± 1.29 


3.54± 1.29 


5.32± 1.23 


T .vn pa 1 


1 


6 


-25.08 ±4.33 


-24.31 ±0.52 


-5.20 ±4.50 


-7.70 ±4. 10 


-8.78 ±4.36 


-3.03 ±4.25 


Lynga 2 


4 


6 


-9.64 ± 1.29 


-9.94 ± 0.60 


-3.88 ± 1.53 


-8.61 ± 1.50 


-8.99 ±1.51 


-2.89 ± 1.52 


NGC 129 


3 




-39.41 ±0.54 




-2.37 ± 1.25 


-1.32± 1.25 


-2.46 ± 1.25 


-1.15 ±1.25 


NGC 381 


2 




-29.76 ± 1.32 




0.63 ± 2.87 


-2. 14 ±2.96 


0.75 ± 2.87 


-2.10 ±2.96 


NGC 457 


2 




-32.99 ± 0.60 




-2.84 ± 1.36 


_0.04± 1.35 


-2.82 ± 1.36 


-0.31 ± 1.35 


NGC 884 


3 




-19.35 ±1.16 




-0.64 ± 1.02 


0.67 ± 1.01 


-0.81 ± 1.02 


0.44 ± 1.01 


NGC 957 


2 




-34.34 ± 1.20 




0.16± 1.35 


-1.15 ± 1.45 


0.54 ± 1.36 


-1.03 ±1.44 


NGC 1513^ 


1 




-15.13 ±1.04 




0.10 ±2.70 


-5.40 ±2.60 


3.60 ±2.66 


-4.02 ± 2.64 


NGC 1528 


5 




-16.90 ±0.48 




1.29 ±0.82 


-0.94 ±0.84 


1.58 ±0.83 


0.14 ±0.83 


NGC 1662 


6 


20 


-12.16± 1.35 


-11. 95 ±0.86 


-1.54 ±0.73 


-2. 17 ±0.73 


1.77 ±0.73 


-1.98±0.73 


NGC 1960 






-17.83 ±0.99 




50±0 51 


_4 5f)±0 55 


4.25 ±0.54 


-1.55±0.52 


NGC 2099 


3 




8.51 ±0.92 




4.53 ± 1.17 


-7.43 ± 1.21 


8.61 ± 1.20 


1.25 ± 1.18 


NGC 2215 


3 


4 


-10.98 ±0.68 


-10.74 ±0.66 


5.07 ± 1.29 


-6.35 ± 1.35 


5.91 ± 1.35 


5.58 ± 1.29 


NGC 2264 


1 


9 


24.06 ±2.91 


24.69 ± 0.98 


-0.40 ± 2.00 


-1.20 ±2.00 


1.16±2.00 


-0.51 ± 2.00 


NGC 2301 


2 




8.25 ±1.14 




-0.55 ±1.13 


-6.20 ± 1.13 


6.20± 1.13 


-0.58 ± 1.13 


NGC 2323 


5 


6 


14.65 ± 1.32 


14.70 ± 1.29 


1.50 ±0.93 


-3.93 ±0.99 


3.82 ±0.99 


1.76 ±0.93 


NGC 2353 


3 


7 


20.16±0.61 


1 8 78 ± 1 44 


-2 75 ± 90 


-2 1 8 ± 89 


2.37 ±0.89 


-2.58 ±0.90 


NGC 2354 


6 


7 


32.89 ±0.39 


32.86 ±0.38 


-9.27 ± 1.25 


-2.13 ± 1.12 


3.74± 1.12 


-8.74 ± 1.25 


NGC 2423 


20 




22.48 ±0.94 




-0.67 ±0.39 


-3.04 ±0.39 


3.08 ±0.39 


-0.43 ±0.39 


NGC 2437 


18 


19 


46.99 ± 1.03 


46.92 ± 1.03 


-4.75 ±0.47 


0.09 ±0.47 


0.27 ± 0.47 


-4.75 ± 0.47 


NGC 2447 


13 




22.38 ± 0.20 




-5.23 ±0.61 


4.33 ±0.58 


-3.69 ±0.58 


-5.70 ±0.61 


NGC 2482 


4 




34 88 + 1 31 




-2 05 ± 92 


4 49 ± 85 


^ 26 ± 85 


-2 49 + 92 


NGC 2516 


5 


57 


23.80 ±0.71 


23.32 ±0.26 


-4.20 ± 0.59 


9.68 ±0.59 


-8.64 ±0.59 


-6.07 ± 0.59 


NGC 2527 


5 


10 


38 73 ± 56 


39 17±0 67 


-3.92 ± 1.40 


6 55 + 1 23 


-6. 14 ± 1.23 


-4.54 ± 1.40 


NGC 2539 


6 


8 


28.36 ±0.41 


28.34 ±0.39 


-3.82 ±0.85 


-3.35 ±0.82 


3.51 ±0.82 


-3.68 ±0.85 


NGC 2546 


7 


9 


18.74 ±0.57 


21.12± 1.32 


-3 .77 ±0.78 


3.74 ±0.72 


-3.31 ±0.72 


^.15 ±0.78 


NGC 2547 


3 


6 


14.22± 1.28 


1 5 65 + 1 26 


_6 74-1-0.73 


3 54-1-0.72 


-2.57 ±0.72 


-7. 17 ±0.73 


NGC 2548 


10 




8.39 ±0.34 




0.61 ±0.48 


1.43 ±0.50 


-1.44 ±0.50 


0.59 ±0.48 


NGC 2567 


3 


6 


35.63 ±0.67 


36.23 ±0.74 


-3.83 ± 1.90 


1.78 ± 1.69 


-1.48 ± 1.69 


-3.96 ± 1.90 


NGC 2579 


2 


5 


1.50 ±0.44 


1.65 ±0.38 


-2.47 ± 1.42 


1.44 ± 1.39 


-1.23 ± 1.39 


-2.58 ± 1.42 


NGC 2669 


1 


7 


20.56 ±0.62 


20.92 ± 0.48 


-3.30 ± 1.30 


6.20 ± 1.30 


-6.10± 1.30 


-3.48 ±1.30 


NGC 2670 


1 


7 


15.74 ±3.44 


17.20±0.83 


-6. 10 ±2. 10 


4.80 ±2.00 


-4.62 ±2.00 


-6.24 ±2. 10 


NGC 2682 


10 


33 


33.67 ±0.42 


33.84 ±0.33 


-7.87 ±0.61 


-5.60 ±0.59 


5.57 ±0.59 


-7.88 ±0.61 


NGC 2925 


2 




-0.25 ± 2.14 




-12.79± 1.55 


3.60 ± 1.45 


-5.38 ± 1.45 


-12.15 ± 1.55 


NGC 3680 


10 


11 


1.06 ±0.36 


1 04±0 35 


_5 75-1-0 61 


1.15±0.57 


-4.45 ±0.59 


-3.81 ±0.60 


NGC 5281 


2 


6 


-19.62 ±0.57 


-18.52 ±0.75 


-9.35 ±1.67 


-0.35 ± 1.59 


-8.25 ± 1.65 


4.40 ± 1.61 


NGC 5316 


8 




-13.55 ±0.48 




-5.27 ± 1.06 


-0.79 ± 1.00 


-4.80 ± 1.04 


2.30 ± 1.02 


NGC 5460 


5 




-8.63 ±2.00 




-5.91 ±0.56 


-2.01 ±0.59 


-5.59 ±0.58 


2.78 ±0.58 


NGC 5617 


3 


6 


-35.95 ±0.80 


-36.60 ±0.86 


1.15 ±2.96 


1.79 ±2.74 


2.11 ±2.83 


0.27 ±2.87 


NGC 5662 


2 




-14.46 ± 2.84 




-2.20 ± 1.72 


-4.77 ± 1.62 


-5.15 ± 1.66 


-1.05 ±1.69 


NGC 5822 


13 




-29.46 ± 0.49 




-8.57 ± 0.95 


-8.85 ±0.88 


-11.69 ±0.89 


3.88 ±0.94 


NGC 5823 


2 


g 


-30.09 ±3.23 


-30.05 ± 0.79 


-4.69 ± 3.24 


-0.60 ±2.98 


-2.57 ± 3.03 


3.97 ±3.19 


NGC 6025 


6 




16.10 ±0.97 




-1.92 ±0.94 


-2.42 ± 0.93 


-2.73 ±0.93 


1.45 ±0.94 


NGC 6031 


2 


4 


-6.23 ±1.76 


-2.48 ± 1.34 


-1.65 ±2.38 


-9.92 ±2.21 


-10.05 ±2.21 


0.13 ±2.38 


NGC 6067 


4 




-39.83 ±0.42 




-2.69 ± 1.10 


-1.76± 1.10 


-2.09 ± 1.10 


2.44 ± 1.10 


NGC 6124 


10 


11 


-19.89 ±0.99 


-19.87 ±0.98 


-0.07 ±0.51 


-3.42 ±0.52 


-3.42 ±0.52 


-0.16±0.51 


NGC 6134 


7 


9 


-27.30 ±0.67 


-27.41 ±0.59 


-0.11 ± 1.12 


-6.97 ± 1.09 


-6.96 ± 1.09 


-0.38± 1.12 


NGC 6167 


3 


5 


-20.53 ±0.73 


-21.32 ±0.41 


-1.94± 1.35 


-1.55 ± 1.31 


-1.63 ± 1.31 


1.86± 1.35 


NGC 6250 


1 


5 


-7.98 ±0.99 


-8.04 ±0.81 


13.50± 1.70 


-11.20 ±1.70 


-1 1.60 ± 1.70 


-13.15± 1.70 


JNLiL. Dzol 


D 


1 1 


— D. j J ± t). /4 


—O.JO it U.oU 


— Z.OD ± v). /j 


—J. DO dz U. /o 


— J. j4 ± \). /O 


3.01 ± 0. /5 


NGC 6405 


6 


10 


-8.27 ±0.45 


-7.02 ±1.53 


-1.49 ±0.63 


-6. 10 ±0.64 


-5.89 ±0.64 


2.18±0.63 


NGC 6416 


6 


10 


-27.02 ± 1.58 


-27.52 ± 0.94 


-0.36 ±0.91 


-0.14 ±0.99 


-0.09 ±0.99 


0.37 ±0.91 


NGC 6603 


2 


4 


21.33±0.96 


21.34 ±0.92 


0.98 ± 1.06 


0.51±1.15 


0.39±1.15 


-1.03 ± 1.06 


NGC 6705 


4 


6 


30.48 ± 1.42 


30.87 ±1.14 


-5.38± 1.19 


-0.35 ±1.25 


-0.05 ±1.25 


5.39± 1.19 


NGC 6811 


7 




6.03 ±0.30 




-5.31 ±0.67 


-8.13 ±0.64 


-9.67 ± 0.65 


0.92 ±0.66 


NGC 6866 


2 




12.18±0.75 




-5.52 ±1.17 


-7.97 ±1.09 


-9.67±1.11 


0.61 ±1.15 


NGC 6885 


2 


6 


-1.60 ±0.99 


-1.50 ±0.87 


-2.90 ±1.41 


-6.05 ±1.34 


-6.70 ±1.35 


0.41 ± 1.40 


NGC 7209 


6 




-20.50 ±0.67 




1.81 ±0.69 


-0.04 ±0.66 


1.44 ±0.68 


-1.09 ±0.67 


NGC 7654^ 


1 




-57.39 ±2.14 




0.40 ±5.00 


0.60 ±5.20 


0.55 ±5.02 


0.46 ±5.18 


Platais 1 


3 




-26.73 ±0.40 




-4.50 ±1.09 


-3.78 ±1.08 


-5.88 ±1.09 


0.16 ±1.08 


Ruprecht 119 


3 


6 


-12.77 ±1.34 


-11.21 ±1.34 


0.92 ±1.46 


-1.25 ±1.44 


-1.18±1.44 


-1.00 ±1.46 


Stock 8 


3 




-18.01 ±1.50 




-1.81 ±1.06 


-4.23 ±1.07 


2.90 ±1.07 


-3.57 ±1.06 


Trumpler 10 


2 


8 


27.62 ±3.10 


32. 17 ±0.76 


-12.17±0.81 


8.05 ±0.78 


-7.81 ±0.78 


-12.33 ±0.81 



Bulk cluster parameters unreliable due to mmbership uncertainty. 
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TABLE 13 

Derived NGC 2682 Motions Compared to Previous 

RliSULl'S 



Parameter 


Units 




^members 


Value 


Ref 


Ha COS S 


(mas yr" 


■') 


10 


-7.87 ±0.61 


This Work 




(mas yr" 


') 


30 


-8.62 ±0.28 


Bias et al. (2001) 




(mas yr"" 


•) 


27 


-8.31 ±0.26 


Kharchenko et al. (2005) 




(mas yr" 


') 


10 


-5.60 ±0.59 


This Work 




(mas yr~ 


') 


30 


-6.00 ±0.28 


Dias et al. (2001) 




(mas yr~ 


') 


27 


^.81 ±0.22 


Kharchenko et al. (2005) 


Vr 


(km s"' 


) 


10 


33.67 ±0.42 


This Work 




(km S-' 


) 


13 


32±9 


Scott et al. (1995) 




(km s"' 


1 


104 


33.5 ±0.5 


Malhicii ct al. (l9S6l 



TABLE 14 
Comparison of Derived Cluster RVs 



Cluster 


Members 
3D 3D+RV 


Bulk Vr,3D 


Bulk Vr,3D+RV 
(kms-') 


Other Vr 

(kms~^) 


# 
stars 


AK,30 


AVr.SD+RV 
(kms-1) 


Reference 


Berkeley 86 


2 


- 


-25.54 ± 2.64 




-22.0 ± 7.0 


6 


+3.5 




Forbes etal. (1992) 


IC 2488 


3 


- 


-1.84 ± 1.52 




-2.6 ±0.1 




-0.8 




Claria et al. (2003) 


IC 4651 


9 


10 


-33.34 ± 0.67 


-33.30 ± 0.67 


-31.0 ±0.2 

-30.8 ± 0.3 


14 

44 


+2.3 

+2.5 


+2.3 
+2.5 


Mermilliod et al. (1995) 

Meibom et al. (2002) 


IC 4756 


7 


13 


-24.44 ± 0.75 


-25.08 ± 0.64 


-25.8 ± 0.6 
— zD.U ± U.z 


13 

1 C 


-1.4 

— U.D 


-0.7 
+U.1 


MermiUiod & Mayor (1990) 
vantova et at. \X3y\}) 


JNUL. 12y 


3 


- 


—5yAL ± U.M 




—Jo. 5 ± U.Z 


L 


+u.y 




Mermillloa et ai. (lyo/) 


IN LiU 4 J / 


2 




— oz.yy iL u.ou 




0< 1 Ion 
— ZD.t it o.U 


A 

4 


+/.y 




T ill iit ol l^ QOQ\ 

i^iu et ai. (.tyoyj 


NGC 884 


3 




-19.35 ± 1.16 




-42.5 ± 2.8 


2 


-23.2 




Liu et al. (1989) 


NGC 957 


2 




-34.21 ± 1.22 




-28.6 ± 13.8 


2 


-5.6 




Liu et al. (1989) 


NGC 1662 


6 


20 


-12.16 ± 1.35 


-11.95 ±0.86 


-13.9 ± 0.5 




-1.7 


-2.0 


Grenier et al. (1999) 


NGC 2099 


3 




+8.51 ±0.92 




+7.7 ± 0.9 


30 


-0.8 




MermiUiod et al. (1996) 


NGC 2264 


1 


9 


+24.06 ± 2.91 


+24.71 ± 0.95 


+24.1 ±8.0 


6 


+0.0 


-0.6 


Liu etal. (1989) 


NGC 2354 


6 


7 


+32.89 ± 0.39 


+32.86 ± 0.38 


+33.4 ± 0.3 




+0.5 


+0.5 


Claria et al. (1999) 


NGC 2447 


13 




+22.38 ± 0.20 




+21.7 ±0.7 


11 


-2.1 




Mermilliod & Mayor (1989) 


NGC 2516 


5 


57 


+23.80 ± 0.71 


+23.32 ± 0.26 


+22.7 ± 0.4 
+24.2 ± 0.2 
+23.8 ± 0.3 
+22.0 ± 0.2 


57 
24 
22 


-1.1 
+0.4 
+0.0 
-1.8 


+0.6 
+0.9 
+0.5 
-1.3 


Robichon et al. (1999) 
Terndrup et al. (2002) 
Jeffries etal. (1998) 
Gonzalez & Lapasset (2001) 


NGC 2539 


6 


8 


+28.36 ± 0.41 


+28.34 ± 0.39 


+29.3 ±0.1 


09 


+0.9 


+1.0 


Mermilliod & Mayor (1989) 


NGC 2546 


7 


9 


+18.74 ± 0.57 


+21.12 ± 1.32 


+16.0 ± 2.0 




-2.7 


-5.1 


Hron (1987) 


NGC 2547 


3 


6 


+14.22 ± 1.28 


+15.65 ± 1.26 


+14.4 ± 1.2 




+0.2 


-1.3 


Robichon et al. (1999) 


NGC 2682 


10 


33 


+33.67 ± 0.42 


+33.84 ± 0.33 


+33.5 ± 0.5 
+33.8 ± 1.3 
+32.0 ± 9.0 


104 
04 
33 


+0.2 
-0.1 
+1.7 


+0.3 
+0.0 
+1.8 


Mathieu et al. (1986) 
Prichet & Glaspey (1991) 
Scott et al. (1995) 


NGC 3680 


10 


11 


+1.06 ± 0.36 


+1.04 ± 0.35 


+0.9 ± 0.2 


6 


-0.2 


-0.2 


Mermilliod et al. (1995) 


NGC 5822 


13 




-29.46 ± 0.49 




-29.0 ± 0.7 




+0.5 




Mermilliod & Mayor (1990) 


NGC 6067 


4 




-39.83 ± 0.42 




-39.9 ± 0.2 


10 


-0.1 




Mermilliod et al. (1987) 


NGC 6134 


7 


9 


-27.30 ± 0.67 


-27.41 ± 0.59 


-26.0 ± 0.2 


14 


+1.3 


+1.4 


Claria & Mermilliod (1992) 


NGC 6705 


4 


6 


+30.48 ± 1.42 


+30.87 ± 1.14 


+34.5 ± 1.4 


29 


+4.0 


+3.6 


Mathieu et al. (1986) 


NGC 6811 


7 




+6.03 ± 0.30 




+7.1 ±0.3 


03 


+1.1 




Mermilliod & Mayor (1990) 


Trumpler 10 


2 


8 


+27.62 ±3.10 


+31.91 ±0.73 


+25.0 ± 3.5 


22 


-2.6 


-6.9 


Robichon et al. (1999) 
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TABLE 15 

Comparison of Derived Mean Proper Motions to those of Dias et al. (2001, 2002a) 







Dias et al. (2001,2002a) 






This Study 






Comparison 










Ma cos S 






/^a cos (5 




^f^OL COS 5 




Dias Memb. 


Cluster 


Stars 


iviciriD. 


(mas yr ) 


(mus yr~ ' ) 


Memb 


(mas yr~ ' ) 


(mas yr~^) 


(mas yr~^) 


(mas yr ) 


Conf. 


kpi 
Kej. 


berkeley (So 


AA 

yy 


50 


-A.\±l.l 


/I C _L 
—4.5 ±2.2 


2 


1 1 1 1 

—3.0 ±1.3 


— 4.0± 1.2 


A 3 

—0.3 


0.1 


2 


22 


Collinaer 


1 A 

ly 


12 


-3.9 ±1.9 


0.5 ±1.9 


1 


/I A _L 1 

— 4.U ± 2. 1 


A 1 1 A 

o.U ±1.9 


A 1 
U.l 


-1.5 


1 


4 


Collinder 25o 


25 


13 


-8.1 ±2.2 


—0.8 lb 2.2 


2 


—0.9 zb 2.6 
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Fig. 1. — Comparison of proper motions AfiacosS and '^fJ-S derived from tlie Hipparcos (Baumgardt et al. 2000), Tyclio-2 (Dias et al. 2001, 2002a), 4M 
(Glushkova et al. 1996) and the new UCAC-2 (Dias et al. 2006) surveys, (a) Hipparcos vs. Tycho-2. (b) 4M vs. Tycho-2. (c) Hipparcos vs. 4M. (d) UCAC-2 
vs. Tycho-2. Error bars are the quadrature combination of the uncertainties in the two surveys. Filled triangles denote clusters with best erors (Ae^^^ and 
Aepj < 2.0 mas yr"'). 
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Fig. 2.— (a) 2MASS color-magnitude diagram (CMD) of all UCAC-2 stars within the 25' radius of M67 with Salasnich et al. (2000) 2MASS isochrone 
overplotted. Red points denote stars brighter than magnitude = 13.0 in the UCAC system (approximately equal to the Cousins R band), that are selected to be 
along the cluster's stellar sequence in the CMD while blue points denote stars fainter than 13.0 that lie the cluster's main sequence, (b) Comparison of proper 
motions /1q and /ig derived for M67 stars from the UCAC-2 survey. The red and blue points denote the same stars as in (a). One can see that by adding the 
fainter UCAC-2 data, and thereby changing which survey the proper motion data ai'e primarily derived from, one can actually change the derived bulk proper 
motion by almost 2 mas yr"' in each direction. The black squai'e denotes the measured bulk proper motion from the Dias et al. (2001) survey. 



Open Clusters as Galactic Disk Tracers I 



23 



o 
a 



isi 



<P 1 

9 

< 8 
ao 7 
o ' 

^500 



\- (b) 



-500 



(a) 



~i 1 r 



"T 1 r 



~i 1 r 



H h 



H ^ h 



▲ ▲ A ■ A A '^J 

H 1 1 V- 



(c) 



A ^ : • • 



J I L. 



J I L 



8R^(kpc) 10 



-2 



-4 



~i 1 1 1 — T 1 1 1 1 1 1 1 1 1 1 r 



(d) 



J I I I -A I I I I I I I I I I L 



8 10 

Xge (kpc) 



12 



A A 



H \ h 



12 



Fig. 3. — Properties of the open cluster sample, (a) Plot of cluster age vs. R^c- One can see that most clusters are less than 200 Myr old (grey triangles) though 
a number of old clusters are present in the sample (black circles). Open Squares denote the clusters CoUinder 258, Lynga 1, NGC 1513, NC3C 6250, and NGC 
7654 (see §6.1 & 6.2.3), which we excluded from our sample because our results for these clusters are uncertain, (b) Same as (a), but showing a smaller age 
range, (c) Distribution of Z^c (height above/below the Galactic plane) versus R^c- (d) Xgc, Igc distribution of clusters in this study. 
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Fig. 4. — Comparison of the measured RVs for Fekel (1999, open boxes) and lAU radial velocity standard stars. The dashed line marks an ideal 1-to-l 
correlation. The dotted line is a linear fit to the data. 
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Fig. 5. — Comparison of the color/temperature of the cluster field stars to the selected RV template, (a) Raw selection of stars into "red", "green", and "blue" 
based on the appearance of the spectra, (b) The distribution of stars cross-coirelated vs. lAU RV standards (red) and Fekel (1999) RV standards (blue). 
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Fig. 6. — Comparison of RVs derived for the "green" spectra (which show both the Ca II triplet and the Paschen series features) using red versus blue cross- 
correlation templates. All stars shown have RV measurement errors in the blue less than 10 km s"\ the black squares have RV measurement errors in the blue 
less than 6 km s"' . The black dashed line shows a perfect 1-to-l correlation and the dotted line is a linear fit to the trends, (a) The March 2002 data shows less 
than a 1 km s~' shift, and therefore no offset was applied between the red and blue data, (b) March 2003 shows a -1 km s~' shift, which was applied to the "blue" 
RVs. (c) August 2003 shows a -10 km s"' shift, which was the correction appUed to the "blue" RVs. (d) September 2003 shows a -6 km s"' shift, which was 
applied to the "blue" RVs. Note: the July 2003 run had no stars where the blue template error was less than 10 km s~' . 
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Fig. 7. — Comparison of our measured RVs versus previously published star-by-star high precision velocity results. All stars are color-coded by the observing 
run on which they were observed as follows: red = March 2002; green = March 2003; cyan = August 2003; and blue = September 2003. (a) AVr vs. the V 
magnitude converted from the Tycho (W) magnitude, (b) AVr vs. the 2MASS Kg magnitude, (c) AVr vs. the (B — V) color converted from the Tycho (VtMt) 
magnitudes, (d) AV,- vs. the (V—K,) color, where V is converted from the Tycho (VV) and (Ks) is the 2MASS magnitude, (e) AV,- vs. our measured V,-. (f) Same 
as (e) except with cluster IC 465 1 removed. 
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Fig. 8. — Our measured kmematical data for NGC 2682 (M67). (a) RV distribution, sliown with 2 km s ' biiming, of all stars with RVs between -100 and 
+100 km s"' measured for NGC 2682. (b) Proper motion distribution of all observed stars with proper motion data in the NGC 2682 field having Tycho-2 data 
with error bars shown. 
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Fig. 9. — Steps in the membership analysis for NGC 2682 (M67) for the 1-D RV distribution (V,)- (a) Kernel-smoothed RV distribution for all stars used in 
the analysis of the data for the cluster NGC 2682. (b) Kernel-smoothed distribution for stars not within the cluster radius (Dias et al. 2002b). (c) Probability 
distribution estimated by [(a)-(b)]/(a). (d) 1-D Gaussian fit to (c). The fit to the "cluster" distribution used to determine the membership probabUity (P^) based 
on the spatially-constrained RV data. 
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Fig. 10. — Steps in the analysis of the 2-D proper motion distribution (Pctnaidy t^s) f'"' NGC 2682 (M67). (a) Kernel-smoothed distribution of all stars used 
in the analysis, (b) Kernel-smoothed distribution for stars not selected to be RV members (see Figure 9; < 0.8). (c) The probability distribution given by 
[(a)-(b)]/(a). (d) The normalized 2-D Gaussian fit to the distribution in panel (c). The resulting fitted "cluster" distribution used to determine P^. 
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Fig. 1 1 . — 2MASS color-magnitude diagram (CMD) for NGC 2682 (M67) for stars inside the cluster radius (Dias et al. 2002b). Crosses (x) denote stars that 
we determined to be non-members. Large circles denote stars selected to be members based on both RV and proper motion criteria. Triangles denote stars that 
have > 70% but which do not have Tycho-2 proper motion data available. 
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Fig. 12. — 2MASS color-magnitude diagram (CMD) for all clusters using stars inside the cluster radius (Dias et al. 2002b). Crosses (x) denote stars with 
proper motion data that we determined to be non-members. Large circles denote stars selected to be members based on both RV and proper motion criteria. 
Triangles denote stars that have > 70% but which do not have Tycho-2 proper motion data available. 
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Fig. 13. — Same as Figure 12. 
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Fig. 14.— 



Same as Figure 12. 
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Fig. 15. — Same as Figure 12. 
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Fig. 16. — Comparison of radial velocities AVr to previous studies (see Table 14). (a) AVr plotted as a function of our measured K-; no obvious systematic 
trend is seen, (b) Histogram of AVr showing that, besides the cases of NGC 457, NGC 884, and NGC 957 which all compared here to the results by Liu et al. 
(1989, see §6.2.2), all of our measurements of the bulk RVs of the clusters are within 5 km s~' of all previous determined cluster measurements, and with the 

peak at AVr = km s"' . 



Open Clusters as Galactic Disk Tracers I 35 



20 



10 



N 
20 



10 



1 I I I I I I I I I I I I I I I I I r 

1 



(a) 



:N6250 Col258 

I 1 1 rrrrrn | 




N1513 

^r^'' I \ \ h 



(b) 



-N6250 Col258 



_] I cnEO L 



r-TT^^ \\\\\V\\\V\\\V\-'' IM 1 



N1513 



J I I L. 



-10 





^f^acos6 (mas yr-') 



~i 1 1 r 



~i 1 1 r 



T 



~i 1 1 r 



10 



~i 1 1 r 



20 



10 



N 
20 



10 



(c) 



H \ \ h 




Lyl 



N6250 

H \ \ h 



(d) 



_l I L. 



Lyl 



N6250 



-10 





(mas yr"') 



10 



Fig. 17. — Comparison of the proper motions A^^tcosiS and ^M<5 derived from our study to those of Dias et al. (2001, 2002a, Table 15). Histograms of Afi 
showing that, besides the cases of Collinder 258, Lynga 1, and NGC 6250, all of our reliable measurements of the bulk RVs of the clusters are within 5 mas yr"' 
of Dias et al. (2001, 2002a) study with the peak at A^a = mas yr"' and Afis = mas yr"' . (a) Apia showing clusters having seven or more 3D members from 
our own analysis (blue histogram), while the green histogram denotes clusters with 3-6 3D members, and red histogram showing those clusters with less than 
three 3D members, (b) as (a) with histograms color-coded by membership from Dias et al., with the blue histgram having > 20 members, green 10-19 members, 
and red < 10 members, (c) A/i^ with same color-coding as (a), (d) A^^ with same color-coding as (b). 



